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Resonance  integral  formulas  as  derived  by  Linderberg  and  Seamans  from  the 
Heisenberg  equation  of  motion  for  the  position  operator  ate  incorporated  into 
Intermediate  Neglect  of  Differential  Overlap  (INDO)  and  Neglect  of  Diatomic 
Differential  Overlap  (NDDO)  semi-empirical  methods.  Two-electron  integrals 
consistent  with  the  NDDO  approximation  are  evaluated  theoretically  yielding  an 
approximate  Hamiltonian  containing  no  pure  parameters.  When  applied  to  a series 
of  test  molecules,  results  indicate  that  bond  lengths  and  angles  are  comparable 
to  those  obtained  with  conventional  methods.  Singlet-singlet  configuration  inter- 
action (Cl)  calculations  yield  spectra  which  tend  to  be  blue-shifted  3,000-5,000 


1NDO/S. 


CHAPTER  1 

INTRODUCTION  AND  BACKGROUND 


Such  diverse  processes  as  photosynthesis,  oxygen  transport,  and  DNA  repli- 
cation, processes  of  central  importance  to  life  as  we  know  it,  are  fundamentally 
molecular  in  nature.  Since  much  of  our  modem  day  understanding  of  atoms  and 
molecules  has  come  from  the  study  of  quantum  mechanics,  it  is  of  great  im- 
portance that  we  be  able  to  study  these  processes  within  a quantum  mechanical 

However,  the  computational  demands  associated  with  the  application  of  quan- 

Indeed,  it  was  not  long  ago  that  ab  initio 1 calculations  on  even  small  molecules 
were  thought  out  of  the  question.  According  to  Levine  [1] 


in  the  50's  there  was  a general  belief  that  meaningful  ab  initio  cal- 
culation of  molecular  properties  for  all  except  very  small  molecules 
was  out  of  the  question.  Quantum-chemistry  books  written  in  this 
period  contain  such  statements  as,  “we  cannot  hope  ever  to  make 
satisfactory  ab  initio  calculations  [for  organic  compounds]"  and  "It 
is  wise  to  renounce  at  the  outset  any  attempt  at  obtaining  precise 


solutions  of  the  Schroedinger  equation  fo 
cated  than  the  hydrogen  molecule  ion." 


The  advent  and  rapid  development  of  the  digital  electronic  computer  and  a 
greatly  increased  knowledge  of  the  problem  itself  have  rendered  such  statements 
inaccurate  predictions  of  the  future.  Today,  the  principal  limitation  regarding  the 

is  the  speed  of  the  computer  and  its  available  storage  capacity.  Nonetheless, 
these  limitations  do  exist  and  are  such  that  ab  Initio  calculations  on  rather  large 
molecules  such  as  chlorophyl,  hemoglobin,  and  DNA  are  either  prohibitively 
expensive  or,  in  the  latter  case,  still  beyond  the  reach  of  current  technology. 

The  importance  of  these  molecules  and  the  desire  on  the  part  of  theoreticians 
to  study  them  go  far  in  motivating  the  need  for  approximate  methods  whereby  the 
fundamental  physics  of  the  ab  initio  methods  are  essentially  preserved  but  further 
approximations  are  invoked  which  greatly  reduce  the  scale  of  the  calculations. 
Semi-empirical  methods  then  go  a step  further  by  replacing  certain  terms  in  the 
formulation  with  values  obtained  cither  directly  or  indirectly  from  experiment. 
It  is  this  empirical  aspect  of  semi-empirical  methods  which  accounts  for  their 

The  importance  of  semi-empirical  methods  is  undisputed.  From  the  topo- 


logical  Hiickcl  ir-elcctron  theory  to  today's  all-v 


:thods,  semi-empirical 


treatments  have  greatly  enhanced  our  insight  into  molecular  electronic  structure. 

Despite  the  great  successes  of  semi-empirical  methods  there  are  still  some  very 
valid  criticisms  leveled  against  them.  One  such  criticism  particularly  relevant  to 
this  study  is  that  current  methods  utilizing  a single  model  are  not  capable  of 
accurately  reproducing  a variety  of  properties  over  a broad  range  of  molecules. 
Typically,  a different  model  Hamiltonian  is  invoked  for  different  properties.  This 
is  somewhat  disturbing  in  that  in  nature,  there  is  only  one  Hamiltonian-thc  exact 
one.  One  would  never  argue  for  example,  that  one  set  of  physical  forces  is  at 
work  in  determining  molecular  geometries  and  another,  different,  set  is  at  work 
determining  molecular  transition  energies.  It  follows  then  that  if  a model  is 
a simplified,  yet  reasonably  correct  representation  of  reality,  this  single  model 
should  suffice.  That  it  does  not,  leads  one  to  believe  that  the  model  is  lacking 
and  needs  improvement. 

The  limitations  of  ab  initio  methods  combined  with  the  past  success  and  future 
potential  of  approximate  methods  is  powerful  motivation  for  further  investigation 
of  semi-empirical  models.  In  this  work  two  possibilities  for  an  improved  semi- 
empirical  model  are  investigated.  The  first  involves  advancing  from  the  1NDO  to 
NDDO  integral  approximation  and  the  second  is  concerned  with  the  formulation  of 
two-center,  one-electron  Fock  matrix  elements-those  terms  primarily  responsible 


li-empiricat  models.  Chapter 


Chapter  three  describes  an  alternative  formulation  of  the  two-center,  one-electron 
Fock  matrix  elements  which  is  subsequently  implemented  at  the  1NDO  and  NDDO 
levels.  The  results  of  geometry  optimizations  and  singlet-singlet  configuration 
mteracuon  calculations  using  the  new  methods  are  presented  in  chapter  four  and 

It  is  our  hope  that  this  work  represents  a step  in  the  direction  of  a single, 
semi-empirical  Hamiltonian  capable  of  predicting  molecular  geometries  as  well 
as  electronic  transition  energies  utilizing  the  fewest  parameters  possible. 

Background 

The  Schrddinger  Equation 

One  of  the  central  problems  of  quantum  chemistry  is  the  solution  of  the 


Hi  = Ei  (1) 

where  H is  the  Hamiltonian  operator  for  a system  of  nuclei  and  electrons,  $ is 
the  wavefunction  describing  the  system,  and  E is  the  energy  of  the  system.  Given 


"--pvS-gis^-EgS 

+ gg^  + ££^ 


and  M nuclei,  respectively.  Since  r,j  = (r,  - R^|,  the  third  term  represents 
the  coulombic  nuclear-electronic  attractions.  Similarly,  r,,  = |r,  - r,-|  and 
Rab  = |R/t  - R-al,  so  the  fourth  and  fifth  terms  represent  coulombic  repulsions 
among  the  electrons  and  coulombic  repulsions  among  the  nuclei,  respectively. 


The  Bom-Oppcnheimer  Approximation 


To  completely  describe  a molecule  we  must  describe  the  motions  of  its  nuclei 
and  its  elections.  If  we  neglect  the  motion  of  translation2,  then  we  have  still  to 
consider  a molecule’s  electronic,  vibrational,  and  rotational  motions,  i.e.  those 
motions  which  give  rise  to  the  internal  energy  of  a molecule.  Although  these 
motions  are  in  actuality  coupled  in  a complex  fashion,  it  would  be  a great 
simplification  if  it  were  possible  to  view  these  motions  as  independent  and 
therefore,  separately  soluble  problems. 

Toward  this  end,  we  will  invoke  the  physical  picture  of  electrons  moving 
among  the  nuclei  so  much  more  rapidly  than  the  nuclei  move  themselves  that  we 


2.  The  translation  of  the 


physical  pic 


is  the  essence  of  the  Bom-Oppenheimcr  approximation  [3]  which  allows  us  lo 
effect  a quasi-separation  of  variables  so  that  the  electronic  coordinates  depend 
parametrically  rather  than  explicitly  upon  the  nuclear  coordinates.  The  form  of 
the  molecular  wavefunction  resulting  from  this  approximation  is 

*(ri,R.)  = *ji.(ri)n(R.)  (3) 

where  S2(Rn)  is  the  wave  function  for  the  now  separate  motions  of  the  nuclei  and 

= £e,»e,  (4) 
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ie  Hamiltonian  for  the  e 


%l  = *«.(*) 


(6) 


is  the  electronic  wavefunction  and  £,/  is  the  electronic  energy.  The  molecular 
Hamiltonian  can  now  be  rewritten  as 


-£^+££^7+ 


nergy  plu 


potential  energy  of  the  nude 


form  a potential  within  which  nuclear  rotations,  vibrations  and  translations  take 
place.  The  minima  of  this  potential  energy  surface  correspond  to  stable  molecular 
geometries.  In  the  present  work  we  are  concerned  only  with  the  electronic 
Hamiltonian  and  therefore  from  this  point  on  we  shall  drop  the  “e/"  subscripts 
and  understand  that  the  electronic  problem  is  implied. 

Even  with  the  considerable  help  afforded  by  the  Bom-Oppenhcimcr  approx- 
imation, the  electronic  Schrbdinger  equation  can  be  solved  exactly  only  for  one- 
electron  systems.  Now,  if  the  universe  were  composed  only  of  molecular  species 
such  as  the  H2*.  He*4*,  and  H20*9  there  would  be  little  left  to  do.  However, 

hold  them  together  and  we  need  some  way  of  dealing  with  all  these  electrons. 
The  Independent  Particle  Approximation 

The  many-electron  Schrddinger  equation  cannot  be  solved  exactly  but  the 
one-electron  problem  is  soluble.  It  seems  reasonable  then  to  wonder  whether 
the  many-electron  problem  might  be  solved  if  it  were  treated,  in  an  approximate 
fashion,  as  many  one-electron  problems.  More  explicitly,  what  we  would  like  to 
do  is  treat  the  W-elcctron  problem  as  though  it  were  N independent  one-electron 
problems.  The  most  obvious  way  of  doing  this,  and 
is  to  simply  ignore  electron-electron  interactions. 


the  crudest  approxin 


Wilhin  such  an  independent  particle  approximation  the  total  Hamiltonian  for 
the  system  is  just  the  sum  of  the  one-electron  Hamiltonians.  Thus,  the  electronic 
Hamiltonian  becomes 


H = h,  + h,  + ...  + hp  (8) 

with 

fc-  -V?  + £ji  (9) 

The  wavefunction  corresponding  to  a Hamiltonian  of  this  form  is  a simple  product 

* = Xl(r,)x3(r,) . . . xw(r«)  (10) 

E = e i + £j  + ...  + e*  (11) 

Electron  Repulsions 

It  is  now  possible  to  make  a correction  (to  the  energy)  for  electron-electron 
repulsions  by  assuming  the  independent  particle  wavefunction  is  correct  and 


using  it  to  calculate  the  cxpcctatic 


9 

F-I>+E“  <■» 

and  gives  rise  to  energy  corrections  of  the  form 

/ 03) 

which  represent  the  mutual  repulsion  of  the  two  charge  distributions  |x.(*)|2 
and  IXiO')l*-  While  diis  in  some  measure  corrects  the  energy  and  can  provide 
a qualitative  assessment  of  the  effect  of  electron-electron  repulsions,  the  wave 
function  is  still  a poor  representation  of  reality  and  gives  poor  quantitative  results- 
especially  as  the  number  of  electrons  increases. 

The  Method  of  Hanrce 

In  1928  Hartree  proposed  the  method  of  the  self-consistent  field  (SCF)  [4], 
This  method  assumes  that  the  potential  experienced  by  a particular  electron  is 
a spherically  symmetric  average  of  the  potential  due  to  the  N-l  other  electrons 
so  that  the  exact  two-electron  potential  of  Eq.  (12)  is  replaced  by  an  effective 
potential  given  by 

v=E^=EEJi(')'  <l4> 

■ ■ jV* 

where  the  Coulomb  operator  has  been  defined  as 


(15) 


10 

which,  when  operating  on  another  charge  distribution  represents  the  classical 
interaction  of  two  charge  distributions.  It  is  now  possible  to  express  the  total 
W-electron  Hamiltonian  as  a sum  of  effective  one-electron  Hamiltonians  of  the 
form 

H?  = hi  + V?*  (16) 

wavefunction 

(Ai + *?")■«(■•)  = ««(i)  (17) 

The  difficulty  with  these  equations  is  clear  from  Eq.  (14);  the  effective 
potential  for  the  i-th  electron  depends  on  the  as  yet  unspecified  N-l  wavefunctions 
describing  the  remaining  electrons.  Hartree  overcame  this  problem  by  employing  a 
series  of  successive  approximations.  He  began  by  assuming  a set  of  wavefunctions 
(orbitals)  which  he  hoped  were  reasonably  accurate  and  proceeded  to  calculate  the 
Hartree  operator  for  the  t'-th  electron  by  evaluating  the  integrals  of  the  effective 
potential.  The  resulting  N equations  are  solved  for  N new  orbitals  which  are  then 
used  to  compute  the  Hartree  operator  for  the  next  iteration.  The  process  is  then 
repeated  until  the  Hartree  operator  of  the  i+1  iteration  is  within  a given  tolerance 
of  the  i-th  iteration  at  which  time  the  potential  is  termed  self-consistent.  The 
resulting  orbitals  are  then  used  to  form  a Hartree  product  wavefunction  for  the 


many-clcctron  system.  As  will  be  shown  in  the  next  section,  Hattree’s  method 
failed  primarily  due  to  the  form  of  the  resulting  product  wavcfunction  which 
ignored  the  physical  requirement  of  antisymmctiy. 

Spin  and  The  Pauli  Principle 

The  concept  of  an  intrinsic  angular  momentum  or  "spin”  as  a fundamental 
property  of  electrons  was  proposed  in  1925  by  Uhlenbeck  and  Goudsmit  [5,6] 
to  account  for  the  fine  structure  of  atomic  spectra-the  doublet  sodium  D line  to 
cite  one  well  known  example.  Based  on  the  results  of  the  Stern  and  Gerlach 
experiment  [7,8],  as  well  as  spectral  evidence,  it  was  further  postulated  that  the 
magnitude  of  this  spin  is  5 unit  leading  to  eigenvalues  for  5,  of  ±£.  Thus, 
electronic  wavefunctions  should  depend  not  only  on  the  three  spatial  coordinates 
but  on  a fourth  spin  coordinate  as  well.  It  is  customary  to  designate  this  spin 
coordinate  u and  to  introduce  two  orthogonal  functions  a(u>)  and  0(uj)  such  that 


Wo)  = (m  = 1 


08a) 


<«!/?)  = </»M  = o 


(18b) 


»V>=«(*+i)l«>  = fl<*> 

**!/?>=*(<•  + 1)1/9)  = ||/J) 


(19b) 


(19a) 


(20a) 
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s.l“>  = |l«) 

».i«  = -\w  (2«» 


It  is  a fundamental  principle  of  quantum  mechanics  that  the  measurement  of 
any  physical  property  which  depends  on  the  coordinates  of  a system  of  identical 
particles  should  not  depend  in  any  way  on  the  arbitrary  labeling  of  the  particles. 
For  two  electrons  described  by  ®(Xt(2)Xi(2))i  there  exists  a probability  density 
given  by  |®(Xi(l)Xj(2))|i  which  must  not  depend  upon  the  permutation  of  the 
two  electronic  coordinates.  Thus. 

l®(Xi(l)xy(2))|J  = I®(X>(1)X.(2))|J  (21) 

*(Xi(l)Xi(2))  = *(x,(l)Xi(2))  (22) 

(the  symmetric  case),  or 

*(Xl(l)X/(2))=  -®(Xi(l)Xf(2))  (23) 

(the  antisymmetric  case).  Experimental  evidence  indicates  that  "there  is  no  place 
for  the  symmetric  solution  in  otdinary  spectrum  theoiy”  [9]  and  it  is  a further 
postulate  of  quantum  mechanics  that  electrons  are  described  by  antisymmetric 
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wavefuncdons.  The  well  known  Pauli  Exclusion  Principle  is  a consequence  of 
this  antisymmetry  requirement,  for  if  an  electron  is  assigned  the  same  set  of 

in  accordance  with  Eqs,  (24). 

#[*i(l)*i(2)]  = -*[x,(l)x.(2)]  (24a) 

«(xi(l)xi(2)]  = 0 (24b) 


Slater  Determinants 

Considering  the  preceding  discussion  it  is  apparent  that  the  wavefunction 
assumed  by  Hartree  for  the  SCF  procedure  is  incorrect.  This  was  pointed  out 
by  Slater  [10]  who  suggested  that  the  assumed  SCF  wavefunction  should  be  an 
antisymmetrized  projection  of  Hartree  products  defined  by 

't"p  = Ail>  (25) 

where  A is  the  antisymmetrizer, 

■"-vK?1-1"  0B 

and  the  summation  of  Eq.  (26)  is  over  the  N\  possible  permutations  of  the 
N particles.  The  exponent  p,  the  parity  of  the  permutation,  is  the  number  of 


ages  required  tc 


into  its  proper  or  normal  order  and  the  operator  A is  Ihe  antisymmetrizer.  Such 

the  indistinguishability  of  the  identical  electrons  since  the  permutations  distribute 
them  among  all  the  occupied  orbitals  resulting  from  the  SCF  calculation. 

The  rows  of  a Slater  determinant  are  labeled  by  electron  coordinates  (space 
and  spin)  and  the  columns  by  spin  orbital  labels.  We  see  then  that  the  interchange 
of  any  two  rows  of  the  determinant  (corresponding  to  the  interchange  of  electronic 
coordinates)  reverses  the  sign  of  the  determinant.  Furthermore,  should  any  two 
rows  of  the  determinant  be  identical  (corresponding  to  two  electrons  with  the 
same  space  and  spin  coordinates),  the  determinant  vanishes-a  manifestation  of 

The  Hartree-Fock  Equations 

Working  independently.  Fock  [11]  and  Slater  [12]  put  Hartree's  intuitive 
method  of  the  self-consistent  field  on  a more  rigorous  footing  by  making  use 
of  the  variation  principle.  The  variation  principle  assures  us  that  the  energy 
obtained  from  an  approximate  wavefunction  represents  an  upper  bound  to  the 
exact  ground  state  energy.  Hence  the  energy  becomes  one  measure  of  the  quality 
of  a wavefunction-the  lower  the  energy  the  better  the  wavefunction.  Therefore, 
by  minimizing  the  energy  of  a Hartree  product  by  varying  the  radial  functions,  it 
is  possible  to  arrive  at  precisely  the  equations  of  Hartree's  method.  Furthermore. 
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since  it  is  desirable  lo  maintain  the  orthonormality  of  the  starting  orbitals,  the 

can  be  carried  out  using  Lagrange's  method  of  undetermined  multipliers. 

When  this  same  procedure  is  carried  out  using  a Slater  determinant  rather 
than  a Hartree  product  as  the  approximation  to  the  ground  state  wavefunction,  the 
famous  Hartree-Fock  equations  are  obtained.  These  are  given  by 

/( 1)X.(1)  = |>s.Xl(l)  (27) 

where  / is  the  Fock  operator,  defined  as 

r N 1 

/0)x.(i)=  + <28> 

J being  the  familiar  coulomb  operator  of  the  Hartree  method  while  K is  termed  the 
exchange  operator  and  gets  its  name  from  the  exchange  of  electronic  coordinates 
which  takes  place  during  its  operation.  The  exchange  terms  in  the  Hartree-Fock 
equations  arise  solely  as  a consequence  of  the  antisymmetric  nature  of  the  many- 

and  exchange  integrals  are  defined  according  to  Eq.  (29)  below. 

•W)x.(l)  = J dnxH  2)^«(2)x.(l)  (29a) 

W)X.(1)  = J ^;(2)i(t(l)x.(2)  (29b) 

Equation  (27)  can  be  brought  into  the  form  of  a standard  eigenvalue  problem  by 
means  of  a unitary  transformation  which  diagonalizes  the 


rix  of  Lagrangia 


Liltipliers  yielding 


/(1)X.(1)  = £.X.(1)  (30) 


The  electronic  energy  within  the  Hartree-Fock  method  is  given  by 

e = (*o|ff|*o>  = £<«>  + E E WW>  - WliO  (31) 

where  |*o)  's  O’6  Slater  determinant  formed  from  the  N occupied  Hartree-Fock 
orbitals  and  represents  the  best  variational  approximation  of  the  ground  state 
wavefunction  using  a single-determinant  representation. 

Molecular  Orbital  Theory 

Historically,  the  Hartree-Fock  method  described  in  the  previous  section  was 
applied  to  atoms  and  the  resulting  differential  equations  were  solved  numerically 

Such  a procedure  for  molecules,  however,  was  impractical  for  all  but  diatomic 
molecules.  Furthermore,  numerical  wavefunctions  do  not  provide  a convenient 
basis  for  molecular  property  calculations  and  other  post-Hartree-Fock  procedures. 

The  difficult  part  of  the  Hartree-Fock  procedure  is  the  calculation  of  the  molec- 
ular orbitals.  Addressing  this  question,  Roothaan  [13]  and 


Hall  [14]  sugge 


of  basis  functions  according 


i = 5>,cri  (32) 

coefficient  for  the  i/-th  atomic  orbital  in  the  i-th  molecular  orbital.  Such  an 
expansion  reduces  the  problem  of  calculating  molecular  orbitals  from  that  of 
solving  the  Hartree-Fock  equation  numerically,  to  that  of  solving  a system  of 
linear  equations  algebraically  for  the  expansion  coefficients  c„,.  Additionally, 
the  Roothaan-Hall  procedure  provides  a definite  analytic  fotm  for  the  molecular 
orbitals  which  can  be  conveniently  used  for  the  calculation  of  expectation  values. 
The  use  of  linear  expansions  in  terms  of  atomic  orbitals  in  the  Hartree-Fock 
self-consistent  field  method  constitutes  Molecular  Orbital-Linear  Combination  of 
Atomic  Orbital  (MO-LCAO)  Theory.  Since  its  introduction,  it  has  been,  and 
continues  to  be,  one  of  the  most  pervasive  and  powerful  tools  for  gaining  insight 

If  we  now  substitute  the  expansion  of  Eq.  (32)  into  Eq.  (30)  we  get 

(33) 

Multiplication  on  the  left  of  Eq.  (33)  by  *,,,  followed  by  integration  over  the 
electronic  coordinates  yields 

^•=Ew 


(34) 


i F and  S are  the  square,  He 
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F„.  = J «„(l)/(l)^(l)dn  (35) 

V = Jm  (36) 

If  there  are  N basis  functions  there  will  be  N molecular  orbitals  and  N equations 
of  the  form  of  Eq.  (34),  all  of  which  arc  expressed  in  the  single  matrix  equation 
FC  = SCE  (37) 

which  is  the  well  known  Roothaan-Hall  equation. 

The  overlap  matrix  in  the  Roothaan  equation  arises  from  the  use  of  a non- 
orthogonal  basis  set.  However,  a new  coefficient  matrix  can  be  defined  according 
to  LBwdin’s  symmetric  orthogonalization  procedure  [15]  as 

C'  = Sic  (38) 

which  in  turn  gives 

C = S_ic'  (39) 

)vc  substitution  for  C in  Eq.  (37)  and  multiplying  on  the  left  by 


Using  the  al 
S~i  yields 


S-lFS"ic'  = S-isS-ic'E 


(40) 


ric  orthogonalization 


Equation  (40)  defines  the  new  matrix  F'  and  effects  a symmet: 
of  the  basis  which  results  in  the  standard  matrix  eigenvalue  problem 

F'C'  = EC'  (41) 

(F1  - E)C'  = 0 (42) 

This  equation  has  a nontrivial  (C'y4  0)  solution  only  when  the  determinant  of 
F'  — E vanishes  which  leads  to  the  so-called  secular  equation? 

|F*  - E|  =0  (43) 

The  determinant  in  the  secular  equation  is  a polynomial  of  degree  N and  therefore 
has  N roots  which  are  the  eigenvalues  of  Eq.  (42).  Diagonalizing  F'  simulta- 
neously gives  its  eigenvalues  and  the  LCAO  coefficients  (since  C'  F'C'  diago- 

F according  to  Eq.  (39). 

Restricted  Closed-Shell  Hartree  Fock  Theory 
Since  the  discussion  of  spin  and  the  Pauli  Principle,  x<  has  referred  to  a 
generalized  molecular  (or  atomic)  spin  orbital-a  function  of  the  three  spatial 

3.  The  secular  equation  gets  its  name  from  a problem  in  celestial  mechanics  which  gives 
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to  introduce  a set  of  restricted  spin  orbitals,  such  that  two  electrons  may  occupy 


(44) 


' Wr)  u(p)  / 

and  a single  determinant  wnvefunction  formed  from  this  set  of  spin  orbitals  such 
that  evety  spatial  orbital  is  doubly  occupied-i.e.,  there  are  N electrons  (N  is  even) 
NI2  of  which  have  a spin  and  NI2  of  which  have  0 spin-is  termed  a closed-shell 
wnvefunction. 

The  Fock  operator  for  a closed-shell  system  is  obtained  by  integrating  over 
the  spin  coordinates  which  yields 


/(l)  = M1)  + 53  l2j*(!)  “ K"*1)] 


(45) 


J.(  1)  = j dr j^(2)i-V.(2)  (46) 

Ka( l)0i(l)=  y *j«(2)i- 0,(2)]  0,(1)  (47) 

The  factor  of  two  in  Eq.  (45)  arises  since  an  electron  in  say  the  6-th  molecular 
orbital  will  give  rise  to  a coulomb  interaction  with  each  of  the  two  electrons  in 
the  a-th  molecular  orbital  and  the  single  exchange  interaction  reflects  the  fact  that 
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[here  is  only  exchange  between  electrons  of  like  spin.  The  energy  of ; 
closed-shell  determinant  is 

E = 2 5>|k|«)  + £ [2(a6|ai)  - (a6|ia)) 

= 2X>~  + E[2^-,.*1 


ricted 


(48) 


: = hu  + Y,  |2di.  - A-.-]  (49) 


The  restricted  open-shell  (16]  and  unrestricted  [17]  Hartree-Fock  methods 
will  not  be  discussed  in  detail  here.  This  is  not  to  downplay  the  importance  and 
utility  of  these  methods.  However,  the  closed-shell  case  demonstrates  tire  essential 
features  of  the  theory  and  provides  a less  complex  framework  for  the  discussion 
of  approximate  semi-empirical  Hamiltonians.  No  loss  in  generality  results. 

The  Closed-Shell  Expression  for  the  Fock  Operator 
We  are  now  in  a position  to  derive  an  expression  for  the  matrix  elements  of 
the  closed-shell  Fock  operator.  These  matrix  elements  are  given  by 

WW)  = + E PWJ.W  - (50) 

where  elements  of  the  Coulomb  matrix  for  a-th  molecular  orbital  are 

(/i|J.|i/)  = (p|  J dr2rpl^-ip.\i>)  = (pall'd) 


(51) 


= (pa|ai/) 


(52) 


for  the  exchange  operator.  Combining  Eq.  (50),  Eq.  (51).  and  Eq.  (52)  gives  the 

/nv  = V + E I2<H-a>  - (53) 

Inserting  the  linear  expansion  of  Eq.  (32)  for  the  a-th  MO  yields 

= V + EE CtotZMurWX)  - (HM1 

= V + 2 E E C>'C'"  - jWm] 

= V + G„„  ' (54) 

The  two-electron  contributions,  G„„,  are  further  refined  according  to 

N/‘‘  r 1 1 

G<- = 2 E E G*-G«  [o»m>  - 5(ha*)J 
= E 2 E G*»G«  [<H^>  - 2<Hte)] 

= E^[wM>-|o«r|A<')]  (55) 

where  P\„  are  the  Fock-Dirac  or  first  order  density  matrix  elements  which  for  an 
orthogonal  set  of  basis  functions  correspond  to  the  charge  (one-center  elements) 
and  bond  Older  (two-center  elements)  matrix  elements  defined  as 

pa„  = 2Eca.c;„ 


(56) 


Thus,  the  final  expression  for  the  elements  of  the  Fock 


ft*  = ft,  + £ ft*  [OH-A)  - JMM] 


(57) 


It  is  the  evaluation  and  manipulation  of  the  two-electron  integrals  contributing  to 
the  Fock  matrix  which  constitute  the  major  difficulty  in  carrying  out  MO-LCAO 
SCF  calculations.  The  number  of  these  integrals  (neglecting  symmetry)  is  given 
by 


Thus  a valence-only  minimal  basis  set  calculation  on  water  (6  basis  functions) 
involves  231  two-electron  integrals.  For  a similar  calculation  on  formaldehyde 
(10  basis  functions),  the  number  is  1,340  and  for  benzene  (30  basis  functions), 
the  number  escalates  to  108,343. 

Since  the  formation  of  the  Fock  matrix  (an  N4  procedure)  and  the  manipulation 
of  the  two-electron  integrals  is  such  a formidable  task  in  molecular  calculations, 
quantum  chemists  have  long  sought  approximations  to  the  full  Hartrec-Fock 
equations  which  reduce  their  scale  while  still  providing  accurate  results  within 
a physically  meaningful  model.  These  approximations  generally  seek  to  reduce 
the  number  of  two-electron  integrals  by  neglecting  a majority  of  them  in  some 
systematic  way.  When  the  neglect  of  some  two-electron  integrals  is  combined 
with  a replacement  of  other  one-  and  two-electron  integrals  by  values  taken  from 
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Since  almost  all  approximate  methods,  including  those  to  be  developed  in  this 
work,  are  semi-empirical  in  nature  we  shall  turn  our  attention  to  a discussion  of 
several  of  these  in  the  next  chapter. 


CHAPTER  2 

SEMI-EMPIRICAL  METHODS 


The  need  for  semi-empirical  methods  can  be  discussed  in  terms  of  the  three 
areas  in  which  it  enhances  the  science  of  quantum  chemistry.  These  are  1)  the 

state-of-the-art  calculations  on  systems  too  large  to  be  attempted  with  ab  initio 
methods;  and  3)  the  ability  to  line  tune  a semi-empirical  method  to  the  point  that 
it  can  predict  experimental  results  more  accurately  than  a strictly  ab  initio  method. 

Even  though  computers  are  growing  ever  faster  and  computing  time  is  be- 
coming correspondingly  less  expensive,  it  cannot  be  said  that  the  need  for  semi- 
empirical  methods  is  just  as  rapidly  diminishing.  It  is  important  that  quantum 
mechanical  methods  should  be  routinely  available  not  only  to  the  theoretical  re- 
searcher but  also  to  the  industrial  chemist  and  the  college  sophomore.  For  this 
to  be  possible,  the  methods  must  be  relatively  small-scale  and  inexpensive.  It 
is  semi-empirical  methods  which,  to  an  appreciable  degree,  meet  this  criterion. 
For  they  reduce  many  interesting  and  practical  applications  to  a scale  that  can  be 
25 
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performed  with  the  computer  facilities  typically  available  to  the  bench  chemist 
or  the  graduate  student. 

A recent  study  of  the  reaction  center  of  photosynthetic  bacteria  [18]  included 
semi-empirical  spectroscopic  calculations  (singlet-singlet  Cl)  which  reproduced 
the  lowest  excitation  to  within  approximately  300  cm*1  of  the  experimentally 
observed  value.  This  calculation  involved  518  atoms  and  about  1,400  (valence- 

numbers  game  which  should  provide  an  estimate  of  the  scale  of  such  a calculation 

Assuming  10  basis  functions  per  atom  or  5,180  basis  functions  implies 
87,249,408,020,000,  or  roughly  87  trillion,  two-electron  integrals.  Storing  these 
integrals  on  disk  as  IEEE  double  precision  floating  point  numbers  will  require 
8 bytes  for  each  integral  plus  an  additional  16  bytes  for  each  set  of  labels  (4 
bytes  times  4 labels)  for  a total  of  24  bytes  per  integral  or  2,094  trillion  bytes  of 
disk  storage.  Now,  using  1.2  Gbyte  magnetic  storage  disks  which  are  becoming 

drives  to  accommodate  the  two-electron  integrals  alone.  If  we  next  assume  that  we 
have  a very  fast  computer  capable  of  say,  100  MFLOPS  (Million  Floating-point 
Operations  Per  Second)  and  that  the  integrals  can  be  pipelined  very  efficiently 
through  the  floating-point  processor,  i.e.  so  that  there  are  no  memory  or  disk 
access  delays,  then  the  SCF  itself  should  take  approximately  700  million  seconds 
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or  about  22  years  per  cycle.  Assuming  we  are  still  interested  in  the  problem, 
we  must  now  transform  the  integrals  in  preparation  for  a Cl  calculation-an  A/5 
procedure  at  best  This  will  take  about  1.1  million  years. 

At  this  point  we  might  begin  to  think  about  the  size  of  the  Cl  matrix  but  1 
think  the  point  is  clcar-if  wc  arc  going  to  apply  quantum  mechanical  methods 
to  systems  such  as  the  photosynthetic  center  mentioned  above,  semi-empirical 
methods  are  the  only  means  possible  given  the  current  level  of  technology. 

Finally,  the  accuracy  with  which  a particular  method  or  model  predicts 
observations  is  not  to  be  overlooked.  For  example,  assuming  the  very  large  ab 
initio  calculation  of  the  above  example  could  be  carried  out,  evidence  indicates  that 
the  result  would  not  be  very  accurate.  This  points  to  another  strong  argument  for 
the  use  of  semi-empirical  methods.  The  ability  to  parameterize  an  approximate 
method  on  experiment  can  lead  to  very  accurate  results  for  a specific  physical 

accurately,  say,  the  spectra  of  third  tow  transition  metals,  may  not  do  a very  good 
job  in  other  areas.  Nonetheless,  the  accuracy  attained  by  some  semi-empirical 
models  makes  them  very  useful,  albeit  specific,  tools  in  the  study  of  molecular 

Whatever  else  wc  ask  of  a theoretical  model,  we  want  first  and  foremost 
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for  it  to  provide  us  with  physically  meaningful  results.  FOr  it  to  do  otherwise 
is  at  best  unsatisfying  and  at  worst,  could  lead  to  scientific  disaster  (or  at  least 
embarrassment).  And,  this  in  feet  was  one  of  the  obstacles  to  success  in  the  early 
development  of  valence-only  semi-empirical  models.  Simple  attempts  to  extend 
PPP-like  theories  to  all-valence  basis  sets  were  unsuccessful  since  they  ignored 
the  physical  requirement  of  rotational  invariance. 

It  was  Pople  [19]  and  separately,  Manne  [20,21],  who  realized  that  the 

lead  to  the  same  results  for  different  bases  related  by  a unitary  or  orthogonal 
transformation.  The  validity  of  this  requirement  can  be  seen  by  considering  the 
non-hybridized  atomic  orbital  bases  y>*  and  y>'*  related  by  a unitary  transformation 
T such  that  ft  — y‘'T.  In  the  transformed  basis  the  Hartree-Fock  equations  can 
be  written  in  matrix  form  as 

F'C'  = eS'C'  (59) 

If  we  now  multiply  from  the  left  of  Eq.  (59)  by  T and  insert  the  unit  matrix 
1 = T'T  we  get 

TF'T'TC'  = cTS'T'TC  (60) 

which  reduces  to 


FTC  = eSTC' 


(61) 
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and  implies  lhat  C = TC'.  Finally,  the  equivalence  of  the  two  resulting  sets  of 


</>•  = v>'C  = ^'TC'  = vp'C  = t/.'1  (62) 

This  invariance  is  assured  in  non-cmpirical  methods  by  the  invariance  of  the 
Fock  operator  to  unitary  transformations. 

However,  when  we  begin  making  approximations  which  reduce  the  number 
of  integrals  in  the  full  Roothaan  LCAO  SCF  equations,  it  may  happen  that  these 
approximations  lead  to  results  which  are  not  invariant  to  transformations  which 

1.  Mix  atomic  orbitals  having  the  same  principal  and  azimuthal  quantum  num- 

2.  Mix  atomic  orbitals  of  differing  principal  or  azimuthal  quantum  numbers  on 
the  same  atomic  centers  (hybrid  atomic  orbitals). 

(symmetry  adapted  orbitals) 

One  of  the  characteristics  common  to  all  of  the  ZDO  models  developed  by 
Pople  and  his  coworkers  is  their  invariance  to  the  first  and  second  types  of 

The  Core  Hamiltonian 

Another  characteristic  commonly  found  in  valence-only  semi-empirical  meth- 
ods is  the  removal  of  inner  shell  electrons  by  assuming  they  contribute  an  effective 
electrostatic  potential  which  screens  the  valence  electrons  from  the  field  of  nu- 


clear charges.  In  the  case  that  the  screening  is  100%  effective  we  can 


valence-only 


£j=_IvJ-£i£.  (63) 

where  Z'c  is  the  atomic  valence  charge  of  nucleus  C.  (Note  that  since  e in  atomic 
units  is  unity,  we  must  explicitly  include  the  minus  sign  in  the  equations.)  Thus, 
the  general  expression  for  the  core  Hamiltonian  matrix  elements  is 


= (fusil's) 

= <wl-|v?-£vch ») 

= (Ml  - jV?M  - YfrMcWa)  (M) 

Vfc  = — (65) 


Before  considering  the  various  ZDO  methods  in  detail  it  is  helpful  to  write 
down  the  general  form  of  the  core  Hamiltonian  within  our  assumption  of  an 
orthogonal  sj>  type  atomic  orbital  basis  set.  As  will  be  the  case  throughout  this 
development  there  are  three  cases  to  be  considered,  namely 
1.  p = i/;  A = B:  One-Center,  Diagonal 

v;  A = B:  One-Center,  Off-Diagonal 
i/;  A $ B:  Two-Center 


2.  p* 

3-  P# 


two-center  parts  as  follov 


A = B: 

BcfAfA  = + VaM  - Y.  MVcM  (66) 


b;a„a  = m - iv? + vaWa)  - Y2  mvoWa)  (67) 

C*A 

ng  to  the  purely  one-center  part  of  Hc  as  U and  noting  that  in  an 
one-center,  off-diagonal  kinetic  energy  integrals  vanish,  as  do  the 
off-diagonal  nuclear  attraction  integrals,  leads  to  the  result 


= vAAAA  - Y,  MVcM 


(68) 


s;AyA=  -Y,MVc\»a)  (69) 
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3.  Utr,  AfiB: 

K*.  = Oul  - - X>x|VcM  (70) 


The  Complete  Neglect  of  Differential  Overlap 
Zero  Differential  Overlap 

The  Complete  Neglect  of  Differential  Overlap  (CNDO)  [19]  is  the  most 
approximate  of  the  methods  to  be  introduced  in  this  chapter.  It  is  implemented 
by  applying  the  zero  differential  overlap  approximation  to  all  products  of  atomic 
orbitals  where  n yl  v.  Thus  the  matrix  of  overlap  integrals  is  simplified 
to  the  unit  matrix.  Furthermore,  this  approximation  restricts  the  types  of  two- 
electron  integrals  which  are  nonzero  to  two;  1)  one-center  Coulomb  integrals,  and 
2)  two-center  Coulomb  inlegrals-all  three  and  four-center,  two-electron  integrals 
are  neglected  (as  is  the  case  in  all  of  the  semi-empirical  methods  considered 
in  this  work).  Further  note  that  there  are  no  exchange  type  integrals  included 
at  the  CNDO  level  of  approximation  (this  is  not  true  of  INDO,  and  NDDO,  to 
be  introduced  in  subsequent  sections,  which  do  not  neglect  one-center  exchange 
integrals).  These  implications  are  summarized  in  Eq.  (71). 

S„t.B  = (71a) 

(71b) 


(IIA^cWb^d)  = 


: from  this  firs!  approximation 


33 

If  we  now  incorporate  the  simplifications  that  result 
into  the  Roolhaan  equations  for  the  closed-shell  case,  the  following  expressions 
are  obtained. 

1.  Case  p = v\  A = B: 

fiM/.-.  = +Y.Y.Y.Y,  P'dXc^aWwd) 

= tf/i  A*A + yz  yz  y^  yz  Proxe^A^c^AOD^co^ 

- ^yzyzyzyz 

= + yz  yz  ^cacOmaci/mac) 

- ^ yz  yz  PwAmpaipapa)  to 

2.  Case  p ? w,  A = B: 

fh“- = p**<>» + yz  yz  yz  yz  p-o^a^cWaod) 

- 1 yz  yz  yz  yz  Pn^A^vn^^ 

= Pf  .'A + yz  yz  yz  yz (^aci^o^cd  wa» 

- $ yz  yz  yz  yz  PnAc^A^oDUA^ADiAc^SiM 

= PfA*»  - \ P?.»a{I‘A‘'a\I‘AI'a) 


(73) 


3.  Case  /r  v; 


A*B: 

*iw»»  = H**"»  + S H zm  p°'>*c(i*axcWb'’d) 

- \ Him  H P»»*e  <ftt>ekl>*») 

= + Z Z Z Z P"!>Xc  [(/‘A*c\‘,aaD)6AB&CD6i->t\* 

- PcD^A^claOVB^AD&Bci^hu 

= 5 P^.^AVB\HA»b)  (74) 


Rotationally  Invariant  Two-Electron  Integrals 

There  is  no  question  that  the  above  expressions  greatly  simplify  the  effort 
required  to  calculate  the  elements  of  the  Fock  matrix.  However,  with  this  sim- 
plification comes  a new  difficulty-namely  we  have  introduced  rotational  variance 

Consider  for  example  the  integral  (saa/t|PiDPi0).  Such  an  integral  is 
nonzero.  However,  if  we  now  rotate  the  local  coordinate  systems  of  atoms  A 
atomic  orbital  transforma- 


by,  say  45°  (counterclockwise),  the  following 


p'b  = 7!(p'°+p"')  <75) 

P»  = ^(Pm-P«) 

This  means  lhai  ihc  integral  (sA«/|pt»p«B)-  which  is  nonzero  in  one  local 
coordinate  system,  becomes 

(s^Sx|p'xop',„)  = ^\(*A*A\Pr,P*B)  + 2(s^s^|p,„p„s)  + (s^.!aIp,oP»£,)]  (76) 

in  the  new  coordinate  system.  Since  integrals  of  the  form 

{•a^a^zb^b)  = i|(»yU,t|PMP,«)  - (*A«^|P«»P.B)I  (77) 

are  to  be  neglected  (i.e.  assumed  to  be  zero)  regardless  of  the  choice  of  coordinate 
system,  it  follows  that 

(*/!«,!  IPjbPm)  = (sM«^IPxbPxii)  (78) 

which  is  consistent  with  the  result 

IpxbPxj,  ) + (a.a.|p„p„)]  (79) 

and  which  becomes  a condition  that  must  be  enforced  if  the  CNDO  method  is 
to  be  rotationally  invariant. 

Therefore,  to  enforce  such  a condition  it  is  assumed  that  the  values  of  the 
two-electron  integrals  are  independent  of  the  specific  atomic  orbitals  and  become 


instead  a property  of  the  atom  A or  the  atom  pair  A-C, 

{/‘A^a\I‘A^a)  = 1AA  (80a) 

{l‘A^c\l‘A^c)  = 7 AC  (80b) 

This  satisfies  the  requirement  of  rotational  invariance  where  the  two-electron 
integrals  are  concerned  and  allows  further  simplification  of  the  Fock  matrix 
element  expressions  to 
1.  Case  » = v ; A = B: 

r„„  - - | + ££  Bu 


2.  Case  /i  v\  A = B: 

f,a>a=bm„a-  \p,a. 


Case  fi  I»;  A ^ B: 


- 2P<‘xi'B7xfl 


(83) 


The  Core  Hamihonia 


Having  discussed  so  far  how  the  ZDO  approximation  is  applied  to  two- 
electron  integrals  within  the  CNDO  method,  it  is  now  necessary  to  examine 
the  effects  of  the  ZDO  approximation  upon  the  one-electron  integrals  and  the 
subsequent  consequences  for  the  form  of  the  core  Hamiltonian.  Beginning 
with  the  expression  of  Eq.  (64)  we  note  first  that  all  three-center  nuclear- 
electronic  attraction  integrals  arc  neglected  and  of  the  remaining  one-  and  two- 
center  integrals,  only  those  contributing  to  the  matrix  elements  along  the  diagonal 
are  retained  according  to  Eq.  (84), 

OmIVbIi'c)  = (ZOO)  (84) 


However,  we  have  once  again  introduced  an  approximation  which  is  not 
rotadonally  invariant  by  neglecting  only  the  off-diagonal,  one-center  nuclear 
attraction  integrals.  This  problem  is  corrected  much  as  it  was  for  the  two-electron 
integrals.  Specifically,  they  are  assumed  to  be  independent  of  the  specific  orbitals 
and  become  instead  a characteristic  of  the  atom  pair  involved.  This  approximation 
is  expressed  as 


MVbM  = {/m|Vb|/m)  = Vab  (Ro 


=)  (85) 


Rewriting  the  Fock  matrix  expressions  in  terms  of  the  latest  assumptions  and 
with  the  core  Hamiltonian  now  written  as 

BMA  = VMA  - Y.  v*c  (*6) 

CM 

yields 

1.  Case  it  = v;  A = B: 

= V„A„A  -Yv*c-  \p^.1aa 

= Ol„u-  JlW-u  + £iW&u 

+ E l-Ktc  + £ ^Wr-icl  (87) 

Ci<4  AC 


2.  Case  p ? i/;  A = B: 

H„A„A  = 0 (88) 

F„A,A  = - | PrAtJrtAA  (89) 

3.  Case  p ^ i/;  A £ B:  This  case  should  really  remain  unchanged  until  the 
next  approximation.  However,  we 


should  i 
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the  development,  all  three-center  nuclear  attraction  integrals  are  neglected. 

= (jxA\  - |vs  - V*  - VB\vB)  (90) 


F„,»  = M - iv!  -VA-  VB\uB)  - \p„miab  (91) 


The  Two-Center  Core  Elements 

Finally,  all  two-center  core  Hamiltonian  matrix  elements  arc  considered  to  be 
empirically  determined  parameters, 

= </mI  - ivJ  -VA-  Vb\vb)  = (92) 

leaving 

F„a„  = P„An  - ^P„.ubiab  (93) 

for  the  two-center  Fock  matrix  elements.  The  choice  of  values  for  these  param- 
eters, the  so-called  resonance  integrals,  is  vety  important  since  they  represent 
bonding  interactions  in  a molecule. 

Equation  Summary 

Before  concluding  this  secdon,  we  will  summarize  the  CNDO  Fock  matrix 
expressions.  In  doing  so  we  introduce  the  following  standard  notation, 

Paa  = J2P^ 

which  is  incorporated  into  the  final  expressions  below. 


(94) 


CNDO  Fock  Matrix  Element  Expressions 


1.  Case  p = i/;  A = B: 

F,apa=  U^  + Paajaa-  \pfArAlAA+Y.[Pccy*c-VAC]  (95) 

C/A 

2.  Case  n £ v\  A = 23: 

= - \P^.7AA  (96) 


3.  Case 


J**!  A + B-. 


(97) 


The  Intermediate  Neglect  of  Differential  Overlap 

Since  the  Intermediate  Neglect  of  Differential  Overlap  (INDO)  [221  and 
the  CNDO  methods  are  so  closely  related,  it  is  perhaps  easiest  to  describe 
INDO  in  terms  of  how  it  differs  from  CNDO.  This  difference  is  found  in  the 
treatment  of  the  one-center,  two-electron  integrals.  Whereas  CNDO  neglects 
which  violate  the  ZDO  approximation,  INDO  is  somewhat 


all  such  integrals 


two-electron  integrals  are  treated  as  in  the  CNDO  method,  retaining  only  those 
representing  Coulombic  interactions  and  setting  their  value  to  a rotationally 

While  it  is  clear  that  the  1NDO  method  includes  more  two-electron  integrals, 
and  is  therefore  less  approximate  than  the  CNDO  method,  it  is  important  to 
understand  how  keeping  these  terms  enhances  our  approximate  description  of 
electronic  structure.  For  a brief  illustration  consider  the  3P,  'D,  and  '5  terms  of 
the  (ls)2(2s)2(2p)J  configuration  of  the  carbon  atom.  These  terms  differ  from 
each  other  in  energy  due  to  the  exchange  integrals  of  the  form 

<P.P*|P»P*>  (98) 

Such  integrals  are  neglected  within  the  CNDO  approximation  thereby  causing  the 
various  terms  to  be  degenerate.  The  INDO  method  retains  such  integrals  and  is 

a more  accurate  electronic  description  of  the  carbon  atom. 

To  arrive  at  the  appropriate  Fock  matrix  expressions  we  simply  proceed  as  in 
the  CNDO  approximation  adding  the  missing  one-center,  two-electron  integrals 


summarized 


INDO  Fock  Matrix  Element  Expressions 


1.  Case  n = i>  ; A = B : 

= UfAVA  + YL  \pcciac  - Vac] 

1 (99) 

+ 5Z  5Z  Pc^a  [{pa^aIwa)  - ^a^aWa^a)] 

2.  Case  ft  yl  v\  A = B: 

F,aua  = LE  P-aAa  [{VARAVA)  - liliA^WAl'A)]  (100) 


3.  Case  ft  + ir,  A + B-. 

F*.»  = /W.  " \p.a*b-IAB  (101) 

The  Neglect  of  Diatomic  Differential  Overlap 

The  Neglect  of  Diatomic  Differential  Overlap  [19]  (NDDO)  is  the  least 
approximate  of  the  methods  discussed  and  is  also  the  easiest  to  describe.  Just 
ns  in  the  CNDO  and  INDO  methods,  the  ZDO  approximation  is  applied  to  the 
overlap  yielding  a unit  overlap  matrix. 


SfABB  = K 


(102) 


As  for  the  two-electron  integrals,  the  ZDO  approximation  is  invoked  such  that 
the  differential  overlap  of  orbitals  on  different  centers  is  zero  so  that  the  general 
integral  is 

(pa^cWbvd)  = • (103) 

So  once  again  all  three  and  four-center  integrals  are  neglected.  As  in  INDO  all 
one-center  integrals  are  kept  as  are  all  two-center  Coulombic  integrals.  Unlike 
INDO,  NDDO  retains  two-center,  dipole-dipole  interactions  of  the  form 

{/‘A^bWa^b)  {/'S4*'!*'?4  M 004) 

These  integrals  are  known  to  be  important  for  a proper  representation  of  lone  pair 
repulsions  in  molecules  such  as  hydrazine  and  hydrogen  peroxide. 

Another  major  distinction  between  NDDO  and  INDO  (or  CNDO)  is  found 
in  the  core  Hamiltonian  matrix  elements.  Since  we  are  no  longer  neglecting 
monatomic  differential  overlap,  the  approximation  of  Eq.  (84)  is  no  longer  valid 
and  all  two-center  nuclear-electronic  attraction  integrals  must  be  retained.  The 
expression  for  the  one-center  core  Hamiltonian  becomes 

B;a„a  = U„^A  - Y,  (/Ml Vc\va)  • 005) 

C/A 

However,  the  two-center  core  Hamiltonian  matrix  elements,  (A  yS  B)  , 

include  three-center  nuclear-electronic  attraction  integrals  of  the  form  (i‘aIVc{>'b) 
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which  consequently  are  neglected  (since  it  would  be  inconsistent  to  keep  these  as 
the  only  three-center  terms)  leaving 

= ("4l  - - VA  - Vail's y (106) 

As  for  CNDO  and  INDO  these  terms  are  usually  treated  as  pure  parameters  in 
accordance  with  Eq.  (92). 

It  is  worth  noting  at  this  point  that  the  largest  three-center  terms  entering  the 
Fock  matrix  element  F/l(,  are  of  the  form 

- y.MvcM + y~!  pxciApa^cWb^c)  (107) 

If  we  now  approximate  the  two-electron  integrals  according  to  (see  Eq.  (169), 
pg.  73) 

{MA^cWb^c)  » = (tXAac\vB‘c)  (108) 

and  the  nuclear-electronic  attraction  integrals  as 

(PX I Veil's)  * Z'c(HA*c\vB>c)  (109) 

then  Eq.  (107)  can  be  written 

-'Z(z'c-pcc){t‘AscWB*c)  010) 

This  suggests  that  the  three-center  terms  might  tend  to  cancel  in  all  but  the  most 
polar  molecules  in  which  cases  these  terms  probably  represent  the  largest  source 


NDDO 
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obtained  from  the  full  expressions.  The  deiails  can  be  found  in  appendix  B and 
the  results  are  summarized  below. 

NDDO  Fock  Matrix  Element  Expressions 


1.  Case  ii  = v,  A = B: 

ftMIU  = 

+ E E E P«cAc  [</mAc|/m<7C)  - {/uPfclftl}] 
+ E E P**x*  - ^(iiMiA\\A«A)\ 


2.  Case  /t  w,  A = B: 
FfA„A  = U„A.A 


+ E EE [(fUAclKd'c)  - (P/t|Vc|l'/l)] 

C+A  Xc  «c 

+ EE [{^aWa"a)  - liWA^AVA)] 


3.  Case  n j=  v\  A ? B 

*!•<.»»  = Pha-b  ~ 5 E E PcbxAvavb^aob) 


(113) 


The  NDDO 


i distinguished  from  its  CNDO  and  INDO 


terpans  in  several  ways.  First,  and  most  obvious,  it  includes  more  of  the  one 
and  two-electron  integrals  necessary  for  a correct  physical  description.  Thus,  it 
is  clearly  a less  drastic  approximation  to  the  full  Roothaan-Hall  equations.  Most 
importantly  perhaps,  NDDO  introduces  a greater  measure  of  directionality  into 
the  Fock  matrix  through  the  two-center,  two-electron  integrals  and  the  nuclear 
attraction  integrals  since  these  ate  not  spherically  averaged  as  they  are  in  the 
CNDO  and  INDO  approximations.  For  CNDO  and  INDO,  direction  is  imparted 
only  through  the  resonance  integrals,  and  then  only  indirectly  in  that  they  are  set 
proportional  to  the  orbital  overlap. 

More  rigorously,  the  question  of  whether  or  not  there  exists  a basis  set  for 
which  the  ZDO  approximation  is  valid  or  almost  valid  has  been  studied  by  several 
investigators.  Roby  (23,24]  has  shown  that,  for  the  two-electron  integrals  at  least, 
the  NDDO  approximation  within  a symmetrically  orthogonalized  minimal  basis 
set  is  equivalent  to  the  Rudenberg  [25]  approximation  for  integrals  in  the  non- 
orthogonal  atomic  orbital  basis.  This  conclusion  seems  to  be  supported  by  the 
calculations  of  Cook  [26]  and  Brown  [27].  No  similar  result  could  be  found  for 
the  CNDO  or  INDO  prescriptions  for  the  treatment  of  these  integrals. 

However,  it  should  be  noted  that  these  considerations  in  no  way  validate  any 
approximate  treatment  of  core  Hamiltonian  integrals.  These  must  be  considered 
separately  and  each  on  its  own  merit.  This  we  shall  have  ample  opportunity  for  in 
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sections  dealing  with 


then  again  when 


Parameterization  Of  Semi-Empirical  Methods 


The  expressions  for  the  elements  of  the  Fock  matrix  within  the  CNDO, 
INDO,  and  NDDO  methods  derived  in  the  preceding  sections  are  essentially 
recipes  for  constructing  the  Fock  matrix.  They  specify  which  integrals,  or 

level  of  approximation.  However,  they  say  nothing  about  how  the  values  of  these 

value  for  the  term  7,1  g representing  all  integrals  of  the  form 


within  the  CNDO  and  INDO  methods.  But  where  should  this  value  come  from? 


it  might  be  a good  idea  to  point  out  what  the  choices  are  and  what  the  motivation 
should  be  when  making  these  choices. 

While  it  is  clear  that  the  objective  of  approximate  methods  is  good  results  at  a 
reduced  computational  expense,  there  remains  the  question  of  how  "good  results" 
are  to  be  judged.  Two  obvious  choices  exist.  The  approximate  results  can  be 
compared,  at  least  in  a qualitative  sense,  with  1)  experiment  or  2)  ab  initio 


(114) 


: this  is  1 
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using  either  a similar  or  more  complete  basis.  It  is  the  opinion  of  the  author  that 
the  former  is  the  preferred  yardstick  with  the  following  explanation. 

There  is  little  question  as  to  the  value  of  comparison  with  ab  initio  results 
when  the  goal  is  to  test  the  validity  of  new  or  different  approximations.  However, 
a method  which  agrees  well  with  wrong  ab  initio  results,  no  matter  how  aesthet- 
ically pleasing  the  method  may  be,  will  soon  be  supplanted  in  importance  by  a 
method  which  predicts  and  agrees  with  experimental  results. 

Assuming  that  the  goal  is  to  obtain  results  which  agree  at  least  qualitatively 
with  experiment,  it  is  necessary  to  consider  how  the  values  which  go  into  the 

many  ways  perhaps  the  easiest,  is  that  the  integrals  could  simply  be  evaluated  over 
the  chosen  basis  set.  However,  given  the  approximate  nature  of  these  methods, 
such  a non-empirica!  approach  is  likely  to  give  results  (a  fortuitous  cancellation  of 
errors  notwithstanding)  which  agree  neither  with  ab  initio  results  nor  experiment. 

Two  other  reasonable  alternatives  remain.  The  first  is  to  take  advantage  of 
our  considerable  knowledge  of  molecules  and  select  values  which  agree  with  their 
corresponding,  experimentally  determined  quantities  (although  some  consider  this 
putting  the  cart  before  the  horse  in  the  sense  that  we  end  up  predicting  results  a 
posteriori).  The  other  and  more  pragmatic  approach  is  simply  to  choose  values 
which  provide  the  desired  results  as  determined,  for  example,  by  a least  squares 


efully  selected  set  of 


In  practice  a combination  of  all  three  of  these  alternatives  are  used.  For 
instance,  it  is  common  to  choose  values  for  the  atomic  parameters  with  ref- 
erence to  experimental  ionization  potentials  and/or  electron  affinities;  to  evaluate 

resonance  integrals  somewhat  arbitrarily  as  a function  of  the  overlap. 

Before  proceeding  with  any  detailed  discussion  of  parameterization  schemes, 
it  is  useful  to  list  those  quantities  which  need  to  be  determined  before  an  SCF 
procedure  can  be  carried  out  These  are  summarized  in  table  1. 


Table  1:  Terms  and  integrals  entering  the  CNDO, 


One-center  atomic 
paramters 

One-center,  two-electron 
Two-center,  two-electron 


VAC  </*aIVc!i'a> 

/'a<'aI®a^a>  <fa‘'aI<'aAa> 

Tab  </iai'bI®a^b> 

ft'At'B  3/<a"b 


One-Center,  Two-Electron  Integrals 


The  atomic  two-electron  integrals  formed  from  Slater-type  orbitals  can  be 
evaluated  by  partitioning  the  integral  into  its  angular  and  radial  parts.  When  this 
is  done  the  integral  is  expressed  in  terms  of  sums  of  radial  integrals,  each  term 
of  which  is  multiplied  by  a particular  numerical  factor  arising  from  the  integral’s 
angular  component.  The  radial  integrals  are  known  as  Slater-Condon  [28]  factors 
and  their  angular  multipliers  are  given  by  the  Clcbsch-Gordon  coefficients  [29]. 

These  factors  are  of  great  utility  in  the  parameterization  of  semi-empirical 
methods  because  they  allow  for  a reduction  in  the  number  of  parameters  per 
atom.  For  example,  the  six  unique  two-electron  integrals  arising  from  a valence- 
only  sj>  basis  set  for  a second  row  atom  can  be  expressed  in  terms  of  three 


Slater-Condon  factors  according  to4 

M«)  = {sp„|ap,)=F4](s,s)  (115a) 

(«vM  = |oW)  (115b) 

<P,Pn  |P,P,}  = f°  + ^(P.P)  (115c) 

(PvP-IPpP*)  = F°  - J:F2(p,p)  (115d) 

(PnPv  I PvPa  > = fS(P>  P)  (Hie) 


The  Slater-Condon  F's  and  G’s  can  be  calculated  given  the  radial  component  of 
the  orbitals.  This  is  typically  done  for  F°.  However,  one  of  the  requirements 

4.  Note  that  the  2s  and  2p  orbitals  are  assumed  to  have  the  same  radial  function. 
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of  a method  which  docs  not  neglect  one-center  exchange  is  that  it  accurately 
reproduce  atomic  term  energy  differences.  This  requirement  is  better  met,  in 
general,  by  using  the  so-called  "experimental”  Slater-Condon  factors  which  are 
generated  from  a least  squares  fit  to  experimental  energy  levels.  Thus,  the 

empirical  methods  and  arc  used  in  conjunction  with  ionization  potentials  and 
electron  affinities  in  the  determination  of  the  atomic  core  parameters  to  be 
discussed  in  the  next  section. 

One-Center  (Atomic)  Core  Parameters 

The  one-center  integrals  of  the  core  Hamiltonian, 

016) 

represent  the  kinetic  and  potential  energy  of  an  electron  in  the  orbital  on 
atom  A,  in  the  electrostatic  field  generated  by  A's  valence  charge.  These  tetms 
could  be  evaluated  explicidy  over  the  atomic  orbital  basis  and  adjusted  for  use  in 
semi-empirical  calculations  via  a core  pseudo-potential  [30,31].  However,  they 
are  almost  always  treated  as  empirical  parameters  with  values  determined  from 
atomic  spectroscopy. 

To  do  this  in  a manner  which  does  not  favor  one  spectroscopic  state  of 
an  atom  over  another  it  is  reasonable  to  use  the  average  energy  of  a given 
electronic  configuration.  This  quantity  is  the  weighted  mean  of  the  energies  of  the 


from  the  configuration.  For  example,  the  electronic 


configuration  ls22s22p2  has  the  average  energy 


E(ls22sJ2p2)  . 


£(‘5)+5£('g)+9£(3P) 
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(117) 


In  general,  the  energy  of  an  atomic  term  can  be  expressed  in  terms  of  the  core 
integrals  U^,  and  the  two-elec  iron  Coulomb  and  exchange  integrals.  As  was 
shown  in  the  previous  section,  the  Coulomb  and  exchange  integrals  can  in  turn 
be  expressed  in  terms  of  Slater-Condon  factors.  Thus,  it  is  possible  to  express  the 
average  energy  of  an  atomic  configuration  in  these  same  terms.  The  expression 
is  obtained  by  summing  the  ITs  for  each  electron  and  then  adding  the  sum,  over 
electron  pairs,  of  the  average  pairwise  two-electron  interaction  energies  [32]  i.e. 

E=£%+OJii«y  (ns) 

For  a valence  2s2p  atomic  configuration  the  only  possible  pairwise  interactions 


which  lead  to 

£(2sm2p")  = mUi.,1 


+ nUlr.lp  + ^(n>  + n)(m  + n ■ 

lG‘ - ln(„  - !)F2 


1)F° 


(119a) 

(119b) 

(119c) 


(120) 


for  the  average  energy  of  the  configuration  having  m s-electrons  and  n p-electrons. 

Using  these  configuration  averaged  energies,  the  ionization  potential  or  elec- 
tron affinity  for  an  s or  p electron  is  given  by 

h,(A)  = E,|.(2*ra- V)  - ®/t(2am2p")  (121a) 

h,(A)  = TSM[2sm2p’-')  - EA( 2sm2p")  (121b) 

Aj,(A)  = £/,(2sm2p")  - EA-  (2sm+l2p")  (121c) 

-42j,(A)  = £x(2jm2p“)  - (2am2p“+l)  (121d) 

which  provide  four  equations  for  the  determination  of  five  atomic  parameters. 
However,  we  have  already  noted  that  F°  is  usually  evaluated  analytically  over 
a basis  of  Slater-type  s orbitals  and  the  remaining  F's  and  C's  are  taken  from 
experiment,  making  it  possible  to  determine  U's  either  from  IP’s  alone  or  from 
an  average  of  IP’s  and  HA’s  as  Pople  did  in  his  CNDO/2  model. 

Two-Center,  One-Electron  Integrals 

The  integral  VAg  is  not  difficult  to  evaluate  analytically  over  a basis  of 
STO's  and  this,  in  fact,  was  done  in  the  original  CNDO/1  parameterization. 
However,  CNDO/1  tended  to  produce  bond  lengths  which  were  too  short  and 
correspondingly  too  high  dissociation  energies.  Pople  traced  this  problem  to  a 
so-called  penetration  effect  described  by  the  term  Zg yAg  — VAg  which  arises  if 
the  two-center  terms  of  Eq.  (95)  are  rewritten  as 


{ZgyAg  - VAg)  - Qgy,\g 


(122) 
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Qb  = Zb-  Pbb  (123) 

describes  the  net  charge  on  center  B.  This  deficiency  in  the  method  was  corrected 
in  the  simplest  possible  way  by  setting  VAB  equal  to  Z'biab  thereby  forcing  the 
penetration  term  to  vanish  and  establishing  the  commonly  used  approximation  for 
the  two-center  nuclear-electronic  attraction  integrals. 

Two-Center,  TVo-Electron  Integrals 

Various  schemes  have  been  used  for  these  integrals  and  several  important 
ones  are  briefly  outlined  below. 

Analytical 

In  the  original  CNDO/1  parameterization  the  two-center  Coulomb  integrals, 
1a  B were  evaluated  analytically  according  to  the  formulae  given  by  Roothaan 
[33). 

Mataga-Nishimoto 


The  Mataga-Nishimoto  formula  [34]  is 


(125) 


and  ~jaa  = I a + Aa  is  Ihc  relation  adopted  by  Pariser  [35]  based  i 


A modified  Mataga-Nishimoto  formula  having  the  form 


(126) 


was  employed  by  Ridley  and  Zemer  [36].  The  value  of  /,  was  chosen  to  reproduce 
the  spectrum  of  benzene  and  has  the  effect  of  raising  the  value  of  the  integrals 
above  those  obtained  from  Eq.  (124). 


Ohno-Klopman 


For  the  various  MINDO  parametcrizal 
Ohno-Klopman  formula  [45] 


[•ftsfl  + 2 ( ) ] 


where  iaa  and  ibb  are  average  one-center  Coulomb  integrals. 
Two-Center,  One-Electron  Matrix  Elements 


These  are  the  so-called  resonance  integrals  usually  obtained  from 

{x  = <7,t}  (128) 


where  Ba  and  Bb  are  determined  empirically  via  a least  squares  fit  to  data  chosen 
for  a specific  model  (spectra  versus  geometry,  for  example). 


Summary 


I have  attempted  to  give  on  overview  of  all-valence  semi-empirical  methods 
by  describing  the  fundamental  expressions  obtained  for  the  Fock  matrix  through  a 
consistent  and  rotationally  invariant  application  of  the  ZDO  approximation.  While 
1 have  very  briefly  touched  upon  the  topic  of  how  such  approximate  methods  are 
parameterized,  I have  by  no  means  given  an  exhaustive  review  of  parameterization 
schemes.  Over  the  years  many  such  schemes  have  been  introduced  and  several 

be  rather  superfluous  to  reproduce  another  such  review,  the  interested  reader  is 
referred  to  the  broad  base  of  existing  literature  [46-51]. 

The  CNDO,  1NDO,  and  NDDO  integral  approximations  and  their  correspond- 
ing SCF  schemes  arc  undoubtedly  well  established.  The  neglect  of  integrals  and 
the  introduction  of  parameters  from  atomic  spectroscopy  seem  very  well  founded 
based  on  the  accuracy  of  the  results  obtained.  The  main  weaknesses  in  these 
methods  then  appear  to  stem  from  one  or  all  of  the  following: 

1.  The  formulation  of  the  resonance  integral 

2.  The  evaluation  of  the  nuclear-electronic  attraction  integrals  as  equivalent  to 
the  corresponding  two-electron  repulsion  integral 

3.  Neglect  of  two-center  core-valence  repulsion 

4.  The  formulation  of  the  two-center  Coulomb  integrals 

5.  The  assumption  that  a minimum  basis  set  with  a parameterized  Hamiltonian 
can  accurately  reproduce  experiment 


Of  these,  the  formulation  of  the  resonance  integrals,  based,  as  they  are,  on 
pure  parameters,  seems  a likely  point  of  attack  in  an  attempt  to  improve  our 
approximate  Hamiltonian.  The  very  simple  form  for  &BAVB  chosen  by  Pople  and 
his  coworkers  in  the  original  CNDO/1  parameterization  is  based  on  a formula 
used  by  Wolfsberg  and  Helmholz  [52]  in  Extended  Hiickel  type  calculations  on 
tetrahedral  ions.  The  appeal  of  such  a formula  rests  in  the  fundamental  notion 
that  the  bonding  parameters  should  be  proportional  to  the  overlap  of  atomic 
orbitals.  This  concept  is  consistent  with  the  chemist’s  intuitive  understanding 
of  the  chemical  bond  as  composed  of  electrons  shared  among  atomic  nuclei. 

In  Extended  Huckel  type  theories  one  can  show  that  the  largest  off-diagonal 
elements  of  the  Hamiltonian  should  be  proportional  to  the  overlap  by  using  the 
Mulliken  integral  approximation  [53] 


In  ZDO  theories  a similar  proportionality  is  suggested  by  the  symmetrical  orthog- 
onalization  of  the  basis  using  a binomial  expansion  of  the  overlap;  i.e„ 


(129) 


X‘  = XS'i 


(130) 


S-t=(l  + A)-l  = l — Ia  + 


(131) 
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However,  in  this  case  the  relationship  to  the  Slater  basis  is  uncertain.  Furthermore, 
the  expansion  of  S“i  is  divergent  since  many  of  the  elements  of  S arc  greater 
than  about  0.25  for  an  all-valence  basis  set 

Other  formulas,  most  of  which  are  modifications  of  the  Wolfsberg-Helmholz 
formula,  have  been  introduced.  These  have  usually  consisted  of  additional 
parameters.  For  example,  Del  Bene  and  Jaffe  [54, 55]  introduced  a multiplicative 

ft,,,.  = Al*<98ff„,s,,,  (132) 

where  I<„„  — 1 and  K„  = 0.585.  Ridley  and  Zemer  [36]  used  the  same 
foimula  but  set  K„„  = 1.266  and  noted  that  such  scaling  can  be  justified  from 
an  examination  of  ab  initio  calculations. 

The  general  form 

ft.*..  = ftutfj?  + i5)H.Rab)  033) 

where  I * is  the  appropriate  valence  state  ionization  potential  for  an  electron  in 
orbital  ft  on  atom  A,  was  chosen  by  Dewar  and  his  coworkers  for  MINDO.  In 
MINDO/1  the  function  /(Rah)  was  set  to 

I(Rab)  = Aab  + (134) 

with  the  parameters  Aab  and  Bah  chosen  for  each  pair  of  atoms  to  properly 
reproduce  molecular  heats  of  formation.  For  MINDO/2  and  MINDO/3,  /(Rah) 


59 

was  given  by  a single  parameter,  C, is,  for  each  atom  pair  with  optimal  values 
determined  from  a least  squares  fit  to  a series  of  reference  molecules-  More  recent 
methods  such  as  MNDO  (56),  AMI  [57],  and  PM3  [58,59]  use  similar  schemes 

It  should  also  be  mentioned  that  much  of  the  success  of  the  schemes  utilizing 
these  formulas  has  been  in  the  area  of  spectroscopy  (Eq.  (132)  and  thermodynamic 
properties  (Eq.  (133)).  Other  molecular  properties,  most  notably  geometries,  are 
often  poorly  predicted.  The  most  notable  example  of  which  is  the  tendency  of 
ZDO  based  methods  to  favor  the  more  highly  bonded  structure  of  two  isomers. 
However,  this  problem  is  not  attributed  to  the  resonance  integral  formulations  so 
much  as  the  neglect  of  so  many  two-electron  repulsions  that  the  remaining  inte- 
grals, primarily  nuclear  attraction  integrals,  are  given  too  much  weight.  Whether 
this  deficiency  can  be  compensated  for  by  an  alternative  resonance  integral  for- 
mulation is  doubted. 

More  relevant  to  the  current  discussion  is  the  indictment  handed  down  by 
de  Bruijn  [60]  charging  that  schemes  which  set  0 proportional  to  the  overlap  are 
correct  only  for  diatomic  molecules  or  nearest-neighbor  interactions  in  polyatomic 
molecules.  De  Bruijn  further  states  that  formulas  for  non-neighboring  atoms  must 
include  parameters  which  take  explicit  account  of  the  geometry  of  the  molecule. 
However,  his  case  was  based  on  a divergent  expansion  of  S~  5 and  his  suggestions 
for  alternative  forms  for  0 would  force  a larger  parameter  space  on  the  method 
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and/or  introduce  complications  which  outweigh  their  benefits. 

Nonetheless,  in  cases  where  the  expansion  is  appropriate,  for  example  ~- 
electron  calculations  on  benzene,  it  is  pointed  out  that  the  sign  of  0 between 
non-ncarcst  neighbors  is  impossible  to  predict  and  may  very  well  be  positive. 
This  conclusion  is  supported  in  part  by  the  calculations  of  Chong  (61]  who  found 

the  prescriptions  for  0 used  in  currendy  popular  ZDO  methods  which  suggest 
only  negative  values  (assuming  a positive  overlap).  These  facts  constitute  the 
motivation  for  an  alternative  form  for  the  resonance  integral  which  shall  be  taken 


CHAPTER  3 

RESONANCE  INTEGRALS 
Introduction 

be  subjected  to  certain  boundary  conditions  derived  from  basic  quantum  mechan- 
ical postulates”,  Linderbcrg  [62]  derived  formulas  for  resonance  integrals  within 
the  Pariser-Pople-Parr  [51,  63-67]  method  by  making  use  of  the  equivalency  of 
the  dipole  length  and  dipole  velocity  operators.  These  formulas  gave  values  which 
were  well  within  the  range  of  those  obtained  from  empirical  parameters  in  con- 
junction with  Mulliken  or  Wolfberg-Helmholz-like  formulas  and  proved  to  be  of 
practical  use  in  a broad  range  of  applications  including  electronic  spectra,  spin 
densities  [68-72],  optical  activity  [73,  74]  and  metal  band  structure  [75], 

Linderbcrg  and  Seamans  [76]  extended  the  method  for  valence-only  s,p 
basis  sets  and  pointed  out  that  the  simple  one-to-one  correspondence  between 
momentum  matrix  elements  and  two-center  core  Hamiltonian  matrix  elements 
disappears  due  to  a coupling  generated  by  the  off-diagonal,  one-center  position 
matrix  elements.  Lipinski  and  coworkers  have  incorporated  these  formulas  into  a 
valence-only  INDO  method  which  has  successfully  reproduced  ground  and  excited 
state  properties  of  a wide  range  of  molecules  [77-82].  Lipinski  and  Lesczczynski 
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subsequently  extended  the  method  further  to  include  d orbitals  [83]  and  Rydberg 
orbitals  [84]. 

In  this  work  the  resonance  integral  formulas  of  Lindcrbcrg  and  Seamans  are 
used  for  the  two-center  core  Hamiltonian  matrix  elements  in  an  NDDO  level 
semi-empirical  method.  General  applicability  to  prediction  of  molecular  spectra 
and  geometries  are  investigated  and  comparisons  with  respect  to  INDO  results 

What  follows  is  the  detailed  derivation  for  an  sj>  basis  set. 

Derivation 

The  Hamiltonian 

The  starting  point  of  this  development  is  the  equation  of  motion  for  the 
position  operator  in  the  Heisenberg  formulation  of  quantum  mechanics  where 
the  time  evolution  of  a physical  system  is  viewed  as  change  in  the  dynamical 
variables  with  time  [85],  the  state  vectors  remaining  fixed. 

(‘35) 

If  H is  the  many-electron  Hamiltonian  which  is  assumed  to  be  of  the  form 

H = h + g (136) 

where  h and  g are  the  one  and  two-electron  contributions  to  the  Hamiltonian, 


[r,ff]  = [r,h]  + [r,s] 


(137) 


Wc  expect  that  at  an  exact  level  of  theory  the  dipole  velocity  will  be  a one-electron 
operator  [86]  and  consequently  any  two-electron  contributions  will  vanish, 

|r,S]  = 0 (138) 

While  Eq.  (138)  certainly  holds  for  the  exact  many-electron  Hamiltonian  and  a 
complete  basis,  it  is  also  true  that  the  use  of  an  approximate  Hamiltonian  or  a 
truncated  basis  or  both  will,  in  all  likelihood,  invalidate  Eq.  (138).  And,  in  fact, 
some  authors  have  raised  objections  to  the  Linderberg-Seamans  formulas  based 
on  this  point,  since  in  practice  H is  usually  the  Fock  operator  and  the  basis  sets 
are  rarely  complete.  Nonetheless,  as  Linderberg  pointed  out,  the  derivations  ate 
only  rigorous  if  the  assumptions,  of  which  Eq.  (138)  is  one,  are  accepted. 

In  light  of  Eq.  (138)  we  now  have 

ift£  = [r,h]  (139) 

■'  WIpM  = (WIMIM  (140) 

As  is  often  the  case  it  is  useful  to  expand  the  commutator  above  via  the  familiar 
resolution  of  the  identity  to  obtain 


i OmIpM  = - <auWa>W«,I«'*>  (141) 


q.  (141)  now  becomes  ihe  focal  point  of  this  chapter  since  it 


ments  of  the  position  and  momentum  operators  to  those  of  the  core  Hamiltonian. 
The  following  sections  will  develop  approximations  which  will  simplify  this  rela- 

Position  Operator  Matrix  Elements 

We  are  now  interested  in  matrix  elements  of  the  position  operator  over  an 
orthogonal  set  of  sjj  type  atomic  orbitals.  Rewriting  r as 

r = RM  + r-R/1  (142) 

OmM***)  = toftwR*  + (palr-Rak/t)]  (143) 

where  the  arises  from  the  orthogonality  of  the  basis  set.  it  should  be  noted 
that,  rigorously,  the  two-center  matrix  elements  of  the  position  operator  do  not 
vanish  within  the  NDDO  approximation  since  r is  not  a spherically  symmetric 
operator.  However,  experience  has  shown  that  such  integrals  have  little  effect  on 
the  calculation  of  properties  such  as  transition  moments  and  polarizabilities  [87) 
so  that  neglecting  them  is  not  unjustified.  Thus  Eq.  (143)  becomes,  in  essence,  a 
working  definition  of  ZDO  in  the  current  context  and  as  such  is  another  assumption 
of  the  development.  Substitution  of  Eq.  (143)  into  the  right-hand  side  of  Eq. 


(141)  gives  (see  appendix  D) 

•Wlpkfl)  = (R-x  - Rfl)(('x|fc|"B) 

+ e,  MS,A.A{xA\h\VB)  + WafaM*)) 
-usI^B.ofPxIAIxfl)  + Sfl))  ] 

+ e,  + (144) 

1 

+ e,  \uA(6,A,AMhWB)  + W^M) 

- ufl(ivs.„(px|A|xfl)  + WofpxIMsfl))  ] 

Two-Center  Momemum  Matrix  Elcmenis 

We  now  wish  lo  determine  the  general  form  of  the  two-center  linear  momen- 
tum matrix.  In  order  to  keep  the  treatment  as  simple  as  possible  we  will  consider 
a diatomic  molecule,  A-B  in  a coordinate  system  where  A is  at  the  origin  and 
B on  the  z axis  at  a distance  R from  A.  If  we  now  consider  the  transformation 
properties  the  momentum  operator  and  the  basis  functions  under  the  C®„  point 
group,  we  find  that  the  s,  and  p,  basis  functions  and  the  operators  p,  transform  as 
the  totally  symmetric  (A/)  irreducible  representation  (irrep)  while  px,  py,  p„  and 
p,  transform  as  the  £/  or  tt  irrep.  Noting  that 


Ei  ® Ei  = E?  + A|  + Ai 


{Coo.} 


(145) 


the  following  generalized  matrix  representation  of  p can  be  verified. 

[ V„e;  V„e,  V„e,  V„e,  1 


* (paIpM  = 


V„e, 


LV„e,  V„e,  V„e, 
n analogous  treatment  of  the  core  Hamiltonian  yields 
[ft,  0 0 ft, 

, , 0 ft.  0 0 

{fuel's}  - n o n 


Comparing  vector  components  among  Eqs.  (144),  (146),  and  (147)  the  following 


V„=  -Rft.  + uxft,  - uflft,  (148a) 

V„  = u^ft.  - uBft,  (148b) 

V„  = -Rft.  + uA0„  - “oft.  048c) 

V„=  uA0..-uB0„  (148d) 

V„  = -Rft„  (148c) 

V„  = uxft.  (1480 

V„  = -Rft,  + u„ft.  - “aft,  d48g) 

V„  = -uBft,  (148h) 

V„=  -Rft,  + uxft,  - uflft.  (148i) 


The  equations  Eq.  (148a)  through  Eq.  (148i)  constitute  an  overdetermined  system 
of  equations  for  the  five  0's  in  terms  of  the  nine  V„„'s.  It  has  been  suggested  that 


ihe  0's  determined  from  a least-squares  fitting  procedure.  Linderberg  felt,  and 
the  current  author  agrees,  that  such  a procedure  is  not  very  aesthetically  pleasing. 
Rather,  we  shall  proceed,  as  Linderberg  did,  to  find  relationships  among  the  del 
parameters  which  will  reduce  their  number. 

Reduction  of  Linear  Momentum  Matrix  Elements 

From  the  symmetry  of  Slater-type  (or  any  central  field  type)  p orbitals  we 


PAp>  > = P*IP*  > 


(149) 


V„,  = V„  = V„  (ISO) 

and  when  combined  with  Eqs.  (148e),  (148f),  and  (148h)  gives 


R0.,=  -u  a9«=  ub0„  (151) 

Furthermore  it  can  be  shown  (though  not  without  some  effort-sec  appendix  D for 
a derivation  of  the  first  of  these)  that 


(152) 


(153) 


which  together  with  Eq.  (150)  reduce  the  number  of  V parameters  to  five.  We 
can  now  make  use  of  Eqs.  (148a)— (148i)  and  Eqs.  (151)-(153)  in  a series  of 
successive  eliminations  to  produce  expressions  containing  only  the  five  betas  and 
the  two  atomic  transition  moments. 

Beginning  with  Eq.  (148i)  and  combining  it  with  Eq.  (153)  and  Eq.  (150) 


- = -S0„  - 2f?/3„ 


+ 2V„  = -R0„-2R0„ 


Substituting  Eq.  (148e)  into  Eq.  (155)  git 


= - 2R0„ 


ch.  after  rearranging,  leaves 


Proceeding  now  with  Eq.  (148c),  we  can  rearrange 
substitutions  suggested  by  Eqs.  (157)  and  (152)  as  follows. 

- If  <W«)  d(Rv„) 


(158) 


Using  Eq.  (148b)  to  replace  V, 


fll  A dR 
d(Rj)„) 
Rl  B dR 

Finally,  differentiation  yields 
and  similarly  for 


r)  d\R(uA0„-uB0..)] 

dR UBl 

d(R0xx)  d(R0„) 


d(uA0„) 

dR 


(161) 


For  0„  we  begin  with  Eq,  (148e)  and  then  use  Eqs.  (150),  (148f),  (160) 
as  follows 


P..= 


-iv„ 


Ru*0-’ 

1 d(uAuB0„) 


R dR 


(162) 


Finally,  /?„„  can  be 
and  (162)  to  arrive  at 
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obtained  using  the  previous  results,  namely  Eqs.  (157) 


(163) 


can  be  completely  determined  from 
i key  clement  of  any  method  based  on 
noted  that  comparison  of  Eq.  (162)  w 

I Am  f 

R dR  ' 


which  arises  from  the  nonzero  z component  of  the  gradient  matrix  along  the 
intemuclear  axis,  suggests  the  choice  of 


S„  = -uaub0„ 


(165) 


It  is  now  possible  to  express  all  the  resonance  integrals  in  terms  of  first  and  second 
derivatives  of  the  rr-x  overlap  as  follows  to  obtain  the  final  expressions  below. 


&.= 


-S„ 


(166a) 

(166b) 

(166c) 


(166d) 


(166e) 
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s orbital  and  the  other  both  s and  p functions  and  the  case  where  each  atom  has 
only  an  s function  are  obtained  in  a similar  fashion  and  are  given  by 
A(a)  - B(s,p): 


S„ 


(167a) 

(167b) 


ft. 


Si 

R 


(168) 


CHAPTER  4 

PROCEDURES  AND  RESULTS 


General  Parameterization 

The  general  expressions  for  Ihe  Fock  matrix  elements  within  the  NDDO 
approximation  have  been  previously  described  as  has  the  general  method  for 
parameterizing  semi-empirical  methods.  Therefore,  the  following  discussion  is 
limited  to  the  details  specific  to  the  1NDO  and  NDDO  implementations  adopted 
in  this  work. 

The  atomic  parameters.  were  the  same  for  both  the  INDO  and  NDDO 

of  an  atomic  configuration  is  the  same  for  each  of  the  two  approximations.  The 
values  chosen  were  based  on  ionization  potentials  and  are  essentially  those  of 
Ridley  and  Zemer  [36].  This  set  of  parameters  has  proven  reliable  over  time  and 
there  seems  to  be  little  evidence  that  further  modifications  would  prove  beneficial. 

For  INDO  the  so-called  neglect  of  penetration  approximation  Vab  = Z'bIab 
was  kept  as  it  provides  an  effective  means  of  compensating  for  the  assumption  of 
an  orthogonal  basis.  The  actual  calculation  of  this  term  using  a non-orthogonal 
basis  tends  to  produce  nuclear  attraction  terms  which  are  too  large.  In  the  NDDO 
case  the  spirit  of  the  neglect  of  penetration  was  used,  but  in  a manner  consistent 
72 
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with  the  NDDO  approximation.  Thus,  these  terms  were  calculated  according  to 
the  approximation  of  Geoppert-Mayer  and  Sklar  [88]  ax 

(169) 

where  Jfl(2)sfl(2)  is  chosen  as  a spherically  averaged  charge  distribution  at  center 
B.  It  is  readily  apparent  that  for  CNDO  and  1NDO  this  reduces  to  and  is  consistent 
with  Vab  = Z'giAB- 

The  INDO  two-electron  integrals  are  taken  from  the  semi-empirical  Slatcr- 
Condon  factors  for  the  one-center  case  and  approximated  as 

(pa»b\i»a>>b)  = (>asb\sasb)  = -TAB  (170) 

for  the  two-center  case.  For  NDDO,  all  two-electron  integrals  are  evaluated  ana- 
lytically over  STO-3G  expansions  [89]  of  the  Slater-type  orbital  basis  functions. 
AlmloFs  VMOL  [90]  integral  package  was  modified  and  used  for  the  evaluation 
of  these  integrals. 

Resonance  Integrals 

Figures  1 through  4 plot  resonance  integral  (/?)  values  of  N2  as  a function  of 
the  intemuclear  separation.  In  each  of  the  four  graphs  the  first  curve  (marked  by  an 
asterisk)  corresponds  to  the  standard  Wolfsberg-Helmholz  (WH)  formula  with  0^ 
= -25.0  eV,  fa  = 1261  and  fir  = 0.585.  The  second  curve  (marked  by  a triangle) 
corresponds  to  the  Linderberg-Seamans  (LS)  formula  for  the  appropriate  0.  Given 
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the  seemingly  different  functional  forms  of  these  equations,  the  similarity  in  the 
resulting  curves  is  rather  remarkable.  In  each  case  the  two  curves  intersect  in 
the  vicinity  of  the  bonding  region  (approximately  1.1  A for  Na)  but  have  notably 
different  slopes  and  curvatures  which  should  lead  to  differing  bond  lengths  and 

formalism  for  resonance  integrals  was  developed  by  Linderberg  (62]  for  the  x- 
only  PPP  method  in  which  only  nearest  neighbor  interactions  are  considered.  In 
light  of  figure  4,  its  success  in  this  area  is  not  surprising.  Also  worthy  of  mention 
is  the  fact  that  given  the  close  correspondence  of  0,x  within  the  two  methods, 
the  x-x  interaction  factor  of  0.58S,  found  to  be  crucial  to  the  accurate  prediction 
of  electronic  spectra,  seems  to  be  inherent  in  the  LS  formulas.  Based  on  these 
figures  we  can  expect  the  LS  formulas  to  yield  resonance  integral  values  in  a 
generally  useful  range.  Whether  these  values  will  lead  to  the  accurate  prediction 
of  molecular  properties  within  an  approximate,  semi-empirical  framework  is  the 
subject  of  subsequent  sections  of  this  chapter. 

Before  proceeding  to  examine  this  question  in  more  detail  I would  like  to 
comment  briefly  on  figure  5 which  shows  clearly  that  the  LS  scheme  produces 
values  for  0 as  a function  of  R close  to  the  conventional  schemes  which  set  0 
proportional  to  the  overlap.  This  graph  further  illustrates  the  need  for  different 

in  addition  to  the  0.585  factor  often  used  for  the  x interaction  when  0 is  set 


result  of  this 
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proportional  to  the  overlap.  The  oscillation  in  the  a-a  curve  is  the 
interaction  crossing  the  axis  (figure  3),  where  the  interaction  itself  is  small.  In 
most  cases  this  fluctuation  is  unimportant  and  occurs  inside  of  the  equilibrium 


P„  VS.  Rab 


Figure  1:  Wolfsbcrg-Helmholz  and  Lmdcibcig-Scamans 
8U  values  for  N2  as  a function  of  liucmuclcar  separation. 


p,„  vs.  R, 


pOTvs.RAB 


Figure  3:  Wolfsbcrg-Hclmholz  and  Lindcrbctg-Scomans 


p„  vs.  R„ 


Figure  4:  Wolfsbcrg-Hclmholz  and  Lindcrbcrg-Scnmans 
;lxx  values  for  N2  as  a function  of  intcmuclcar  separation. 


Pnv/Spv  VS.  Ra 


INDQ  vs.  NDDO  Matrix  Elen 


One  of  the  most  important  aspects  in  the  development  of  an  approximate, 

Ionian.  It  has  already  been  mentioned  that  the  exact  calculation  of  nuclear- 
electronic  attraction  integrals  in  the  original  CNDO/1  method  led  to  an  approx- 
imate Hamiltonian  in  which  this  balance  was  not  quite  preserved.  The  conse- 
quences were,  in  general,  the  prediction  of  bond  distances  which  were  too  short 
and,  more  specifically,  to  such  phenomena  as  a bound  first  triplet  state  of  H2. 

Pople  attributed  this  problem  to  the  penetration  of  the  core  of  atom  B by 
electrons  in  orbitals  on  atom  A.  This  can  be  seen  by  rewriting 

£ \pcciAC  - Vac ) 071) 

53  (T’cc  - zc)~iac  + 53  Izc7ac  - VAc]  (172) 

We  then  have 

Pcc  -Zc=  ~(Zc  - Pcc)  = -Qc  (173) 

where  Qc  is  the  net  charge  on  center  C.  Following  Geoppcrt-Mayer  and  Sklar  [88] 
Zciac<  can  be  considered  an  approximation  to  the  nuclear  attraction  integral  Vac- 
In  this  case  the  terms  in  the  second  sum  of  Eq.  (172)  can  be  thought  of  as  an  excess 
nuclear-electronic  attraction,  or  penetration  effect,  which  can  be  conveniently 
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neglected  by  setting  Vac  = Z ciac • Viewed  another  way,  the  assumption  of 
an  orthogonal  basis  is  not  compatible  with  the  analytical  calculation  of  nuclear- 
electronic  attraction  integrals  over  a non-orthogonal  STO  basis.  Calculation  of 
the  nuclear  attraction  terms  followed  by  a correction  for  the  core  pseudo-potential 
and  symmetric  orthogonalization  does,  in  fact,  lead  to  Vac  53  Zciac  [91]  which 
seems  to  adequately  compensate  for  this  error,  a posteriori . 

In  light  of  this  discussion  it  is  reasonable  to  ask  whether  the  analogous 

balance  of  the  semi-empirical  Hamiltonian.  Table  2 gives  a term  by  term 
breakdown  of  the  one-center  Fock  matrix  elements  of  N2  within  the  INDO  and 
NDDO  methods5.  Neglecting  momentarily  the  s-p,  element,  the  spherical  nature 
of  the  INDO  gammas  is  readily  apparent  as  the  Zb~iab  INDO  terms,  which  are  the 
average  of  the  corresponding  NDDO  terms.  It  is  also  apparent  that  the  inclusion 
of  NDDO  type  integrals  provides  a more  accurate  representation  of  the  anisotropy 
of  the  nuclear  attraction  and  two-electron  contributions  to  the  Fock  matrix. 

The  NDDO  s-p:  Fock  clement  which  has  no  one-electron  contribution  at  the 
INDO  level,  is  roughly  four  times  smaller  than  its  INDO  counterpart.  Although 
the  NDDO  integrals  give  rise  to  a large  two-electron  contribution,  this  effect  is 
more  than  compensated  for  by  the  one-electron  Zg(pu\ss)  term  resulting  in  the 
overall  smaller  matrix  element.  These  off-diagonal  elements  however,  are  small 
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compared  lo  the  diagonal  elements  and  their  effect  on  the  resulting  wavefunction 
does  not  appear  to  be  very  important. 

Table  2:  Contributions  to  and  values  of  the  1NDO 
and  NDDO  onc-ccntcr  Fock  matrix  elements  of  N2 

1NDO 

Element  F U Zbtab  C 

s-s  -0.879102  -0.940096  -4.922956  4.983950 

z-z  -0.215190  -0.486953  -4.922956  5.194719 


x-x.y-y  -0,136620  -0.486953  -4.922956  5.273289 


NDDO 

Element F U Zb  Onto) 

s-s  -1.030826  -0.940096  -4.906097 

z-z  -0.387308  -0.486953  -5.114161 

s-z  -0.018901  0.000000  -0.600417 

x-x,y-y  -0,294039  -0,486953  -4.827421 


G 

4.815367 

5.213806 

0.581516 

5.020355 


As  to  the  balance  in  the  Hamiltonian,  the  nuclear  attraction  integral  contri- 
butions to  the  one-center  Fock  matrix  are  comparable,  as  one  would  expect  for 
the  two  approximations.  Relative  to  the  averaged  INDO  values,  the  NDDO  z-z 
element  is  slightly  larger  and  the  x-x  and  y-y  elements,  slightly  smaller.  However, 
the  two-electron  contributions  have  increased  and  decreased  respectively,  which 

in  larger  molecules  there  is  cumulative  effect  due  to  the  sum  over  other  centers 


unbalanced  contribution.  However, 


which  could  lead  to  an 

conclusion  is  supported  by  the  overall  reasonable  geometries  obtained  in  the  case 
of  larger  molecules  in  both  methods. 


Ibblc  3:  Contributions  to  and  values  of  the  INDO 
and  NDDO  two-center  Fock  matrix  elements  for  N2 


-0.469041  -0.410835  +0.058206 
+0.461225  +0.295723  -0.165502 
-0.501348  -0.403199  +0.098149 
-0.479978 -0.256270 +0.223708 


NDDO 

Element Fj^ faw G,,,, 

ss  -0.587544  -0.410835  +0.176709 

oo  +0.564756  +0.295723  -0.269033 

so  -0.609809  -0.403199  +0.206610 

rnr -0.488759 -0.256270 +0.232489 


As  can  be  seen  from  table  3,  the  INDO  one-electron  contribution  to  F„„  has 
been  maintained  (as  givien  by  the  Wolfsberg-Helmholz  formula)  in  our  NDDO 
model.  Thus,  the  differences  are  entirely  due  to  the  inclusion  of  the  NDDO 
integrals.  In  general  the  magnitude  of  each  NDDO  element  is  slightly  greater 


than  its  corresponding  INDO 


In  this  section  are  presented  the  results  and  some  discussion  of  four  different 
groups  of  calculations.  Each  group  is  comprised  of  a common  set  of  molecules 
(chosen  somewhat  arbitrarily)  but  differ  in  the  method  used.  Those  labeled  NDDO 
and  NDDO/M  were  both  carried  out  using  the  implementation  of  NDDO  described 
previously  but  differ  in  the  form  chosen  for  the  resonance  integrals;  the  former 
using  the  WH  type  formula  of  Eq.  (128)  and  the  latter  using  the  formulas  due  to 
Linderberg  and  Seamans  discussed  in  the  previous  chapter.  Corresponding  sets 

accordingly.  Experimental  results,  where  available,  are  also  presented. 

Diatomics 

The  NDDO  bond  lengths  for  the  diatomic  molecules  listed  in  Table  4 are, 
with  the  exception  of  U2,  slightly  shorter  than  their  INDO  counterparts.  Lij  is 
slightly  longer  than  its  INDO  counterpart  but  both  the  INDO  and  NDDO  results 
are  far  enough  from  the  observed  value  that  the  difference  is  hardly  significant. 

The  NDDO/M  and  INDO/M  results  indicate  that  for  homonuclear  diatomics 
these  methods  will  give  generally  shorter  bond  lengths  compared  to  those  obtained 
from  NDDO  and  INDO.  For  the  heteronuclear  diatomics,  CO  and  HF  are  slightly 
longer  while  LiH  is  slightly  shorter. 
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Tabic  4:  Calculated  and  observed  equilibrium 
bond  lengths  (A)  for  some  diatomic  molecules 


Molecule  NDDO  NDDO/M  INDO 


H2  0.697846 

U2  2.16S0 

N2  1.1285 

02  1.1253 

F2  1.1027 

LiH  1.5313 

HF  0.9693 

CO  1.1778 


0.644070  0.697646 

1.8621  2.1490 

1.0801  1.1488 

1.0689  1.1421 

1.0697  1.1140 

1.4927  1.5314 

0.9729  0.9781 

1.2022  1.1953 


INDO/M  Obs.  [92] 
0.643852  0.74172 

1.8377  2.6725 

1.1035  1.094 

1.0908  1.20739 

1.0877  1.435 

1.4887  1.5953 

0.9854  0.9171 

1.2461  1.1281 


RMS  0.1100  0.1618  0.1116  0.1650  N/A 

AVG  0.0279  0.0437  0.0261  0.0409  N/A 


Polyatomics 


Comparing  the  INDO  and  NDDO  results  for  water,  we  see  the  now  expected 
shortening  (slight)  of  the  O-H  bond.  An  analysis  of  the  Fock  matrix  elements 
similar  to  that  discussed  for  N2  indicates  that  this  is  more  due  to  an  increased 
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nuclear-attraction  term  (residual  penetration  effect?)  along  the  bonding  axis  than 
an  effect  of  the  NDDO  two-electron  integrals. 

More  interesting  is  the  change  in  equilibrium  bond  angle  from  102.6°  to 
1 12.13°.  This  I believe  is  a consequence  of  the  NDDO  type  two-electron  integrals. 
Of  the  coulombic  integrals  (an«ol»H«o).  (*HP«ol4HPx0)  and  (shP«J*I1Pk)},  the 
smallest  is  the  first  which  is  the  INDO  ~iab-  Therefore  it  is  reasonable  to  conclude 
that  relative  to  NDDO,  INDO  will  tend  to  underestimate  the  angular  coulombic 
repulsions  leading  to  generally  smaller  bond  angles. 

Table  5:  Calculated  and  observed  equilibrium  bond 
lengths  (A-B,  A)  and  angles  (A-B-C,  Degrees)  for  water 

Bond/Angle  NDDO  NDDO/M  INDO  INDO/M  Obs.  [93] 
O-H  1.0018  0.9488  1.0193  0.9739  0.9571 

H-O-H  112.13  130.7170  102.64  116.15  104.52 


If  we  now  look  at  the  NDDO/M  case  we  sec  that  while  the  bond  length  is  not 
unreasonable,  the  bond  angle  has  increased  to  a point  that  is  no  longer  acceptable. 
The  cause  now  is  a combined  effect  of  the  increased  angular  repulsions  discussed 
above  as  well  as  a tendency  for  the  LS  formulas  to  yield  a too  large  s-s  interaction. 
This  results  in  a lowest  a bonding  with  very  little  p character.  By  scaling  by 
a factor  of  0.75,  additional  p character  is  introduced  into  this  molecular  orbital 
reducing  the  bond  angle  (for  the  NDDO/M  case)  to  a more  acceptable  106°  with 
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a bond  length  of  1.102  A.  Scaling  of  the  s-a  interaction  by  1.3  introduces  still 

bond  length  of  0.974  A and  bond  angle  of  102.5°. 

HCN 

The  geometry  of  HCN  is  well  reproduced.  This  is  not  easy  to  understand 
given  that  the  equilibrium  bond  length  for  CO  is  calculated  to  large  while  that 
of  N2  is  calculated  to  small. 

Table  6:  Calculated  and  observed  equilibrium  bond  lengths  for  HCN 
Bond  NDDO  NDDO/M  INDO  INDO/M  Obs.  (931 
C-H  1.0619  1.0101  1.0752  1.0243  1.0632 

C-N  1.1666  1.1273  1.1882  1.1493  1.1552 


Tables  7 and  8 give  the  optimized  bond  length  and  angle  (in  the  C2v  case) 
for  the  C2V  and  D31,  geometries  of  ozone.  For  C%,  ozone,  as  for  water,  there 
is  a slight  increase  in  equilibrium  bond  angle  accompanied  by  a slight  decrease 
in  bond  length  on  going  from  INDO  to  NDDO.  This  is  consistent  with  previous 
observations  and  a similar  analysis  can  be  applied. 


lengths  (A-B) 


nglcs  (A-B-O 


for  ozone  (Cj.) 

Bond/Angle  NDDO  NDDO/M  INDO  1NDO/M  Obs.  [94] 
0-0  1.1643  1.1057  1.1747  1.1206  1.278 

0-0-0  122.46  136.24  120.34  13Z79  116.49 


The  LS  resonance  integrals  give  an  unacceptably  large  bond  angle  for  both 
the  INDO/M  and  NDDO/M  cases.  Unlike  H2O  however,  simply  scaling  the  s-s 
interactions  leads  to  only  a moderate  improvement  reducing  the  angle  by  about 
6°  and  increasing  the  bond  length  by  0.02  A.  Additional  scaling  of  the  s-a 
interactions,  also  by  0.75,  and  the  0-0  interactions  by  1.125  gives  a bond  angle  of 

used  to  improve  the  H20  results  and  indicates  that  INDO/M  and  NDDO/M  results 

integrals. 

Table  8:  Calculated  equilibrium  bond  length  for  ozone  (D31,) 

Bond  NDDO  NDDO/M  INDO  INDO/M 

0-0  1.2415  1.1941  1.2470  1.2122 


An  important  aspect  of  INDO  calculations  on  ozone  is  the  extent  to  which  they 


incorrectly  favor 


1-ring  E>3h  geometry 
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C2«  geometry.  Table  9 gives  the  energies  of  the  optimized  Djh  and  C2v  structures 

between  the  4tr-  and  6!r-electron  states  that  comprise  the  C2,  and  D3h  structures 
lowers  the  C2v  state  by  0.118  a.u.  Thus  the  1NDO/M  and  NDDO/M  procedures 
arc  capable  of  correctly  predicting  the  C2,  form  lower  in  energy  where  the 
conventional  INDO  and  NDDO  methods  are  not. 


ozone  within  the  NDDO.  NDDO/M,  INDO.  INDO/M  methods. 
Geometry  NDDO  NDDO/M  INDO  INDO/M 

D3h  -51.241711  -50.472626  -50.455100  -49.636680 

C2,  -51.051192  -50.440249  -50.222322  -49.572490 

Difference  0.190519  0.032377  0.232778  0.064190 


Formaldehyde  (H2CO) 


A generally  observed  feature  of  geometries  obtained  from  ZDO  methods  is 
that  the  results  tend  to  become  more  accurate  as  the  size  of  the  molecules  increase. 
This  effect  is  usually  attributed  to  a fortuitous  cancellation  of  errors  as  more 

the  formaldehyde  results  this  trend  begins  to  appear  as  the  results  for  all  four 
methods  are  quite  reasonable. 


Tabic  10:  Calculated  and  observed  formaldehyde 
equilibrium  bond  lengths  (A-B)  and  angles  (A-B-C) 

Bond/Anglc  NDDO  NDDO/M  INDO  INDO/M  Obs.  [95] 
C-O  1.2427  1.1963  1.2604  1.2140  1.2077 

C-H  1.0808  1.0353  1.0945  1.0504  1.1159 

H-C-H  118.29  113.54  118.23  112.15  116.50 

H-C-O  120.86  123.23  120.88  123.92  121.75 


Acetylene,  Ethylene  and  Ethane 


Tables  1 1 , 1 2,  and  1 3 give  the  calculated  and  observed  geometries  of  acetylene 
(C2H2),  ethylene  (C2H4),  and  ethane  (C2Hj),  respectively.  These  molecules 
provide  familiar  examples  of  saturated  and  unsaturated  carbon -car  bon  bonds. 
Semi-empirical  methods  are  often  parameterized  to  match  these  experimental 
distances  as  closely  as  possible.  The  experimental  range  of  bond  distances  in 
going  from  the  carbon-carbon  triple  bond  of  acetylene  to  the  carbon-carbon  single 
bond  of  ethane  is  0.3390  A.  This  range  is  not  well  reproduced  by  any  Hartree- 
Fock  method,  the  problem  being  that  the  C-C  bond  length  in  ethane  is  predicted 
to  be  too  short. 

INDO/M  gives  a range  for  these  bonds  of  0.2799  or  slightly  longer  than 
the  0.2774  predicted  by  INDO.  The  NDDO  methods  extend  the  calculated  range 


0.01  A to  0.2823  for  NDDO/M  and  0.2851  for  NDDO.  While 


this  is  only  about  1%  improvement,  it  is  at  least  movement  in  the  right  direction. 
Overall,  agreement  among  each  of  the  methods  with  experiment  is  generally  good. 
Comparing  the  methods  among  themselves  show  that  the  trend  of  NDDO/M  < 
INDO/M  < NDDO  < INDO  maintained  for  bond  lengths  and  that  all  four  methods 
reproduce  the  bond  angles  in  these  molecules  rather  well  (all  within  less  than  4° 
of  experiment).  Each  of  the  methods  correcdy  reproduces  the  increase  in  the  C-H 
bond  length  from  acetylene  to  ethane. 


Bond  NDDO  NDDO/M  INDO  INDO/M  Obs.  [93) 
C-C  1.1847  1.1683  1.2024  1.1904  1.2030 

C-H  1.0609  1.0074  1.0751  1.0198  1.0598 


lengths  (A-B)  and  angles  (A-B-C)  for  ethylene 
lond/Angle  NDDO  NDDO/M  INDO  INDO/M  Obs.  [931 
C-C  1.3116  1.3063  1.3233  1.3107  1.3390 

C-H  1.0776  1.0263  1.0925  1.0452  1.0859 

C-C-H  123.06  122.76  123.4  123.52  121.19 

H-C-H  113.87  114.47  113.21  112.95  117.63 
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Tabic  13:  Calculated  and  observed  equilibrium 
bond  lengths  (A-B)  and  angles  (A-B-Q  for  ethane 

Bond/Angle  NDDO  NDDO/M  INDO  INDO/M  Obs.  [96] 
C-C  1.4698  1.4506  1.4798  1.4703  1.5429 

C-H  1.0846  1.0415  1.0985  1.0570  1.1019 

C-C-H  111.09  110.80  111.49  111.33  109.64 

H-C-H  107.80  108.11  107.38  107.55  109.30 


Although  experimental  evidence  for  the  existence  of  borirane  has  been  found 
[97-100],  the  molecule  itself  has  not  yet  been  isolated  and  characterized.  Exten- 

form  of  C2H5B  labeled  (I)  in  Figure  6 is  0.0868  a.u.  lower  in  energy  than  the 
planar  structure  (II)  of  the  same  figure. 


(i)  <n) 
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Semi-empirical  INDO  mclhods  are  notorious  for  overestimating  the  stability 
of  small,  closed-ring  structures  such  as  (I)  relative  to  a corresponding  open-ring 
structure  (II).  This  is  also  the  case  with  borirane.  INDO  calculations  lead  to  the 
conclusion  that  the  closed-ring  structure  (I)  is  0.21  a-u.  more  stable  than  the 
open-ring  structure  (1). 

In  order  to  determine  whether  or  not  NDDO  or  the  use  of  LS  resonance 
integrals  (or  a combination  thereof)  is  capable  of  correcting  this  error  in  part  or 
in  whole,  the  closed  and  ring  structures  (I  and  II)  of  borirane  where  optimized 
using  these  methods.  The  results  are  presented  in  Table  14  and  indicate  that  the 
tendency  to  favor  small-ring  formation  is  still  quite  strong.  Furthermore,  NDDO 
alone  seems  to  intensify  the  problem  rather  than  improve  it  relative  to  INDO.  This 
is  somewhat  surprising  since  it  might  be  expected  that  the  anisotropic  character  of 
the  NDDO  two-electron  integrals  would  tend  to  redistribute  electronic  repulsions 

this  is  not  observed. 


Table  14:  Calculated  borirane  energies  (a.u.)  for 
both  the  closed  (1)  and  open  (11)  structures 

Structure  NDDO  NDDO/M  INDO  INDO/M 

I -20.8121 18  -18.521706  -19.846130  -17.539277 

II  -20.550769  -18.362983  -19.636571  -17.427151 

Difference  -0.261349  -0.158723  -0.209559  -0,112126 
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There  is  some  cause  for  optimism  however,  in  lhat  similar  calculations 
utilizing  Ihe  LS  resonance  integrals  reduce  tile  extent  to  which  the  closed  ring  form 
is  favored  by  0.097  a.u.  for  the  INDO  case  and  0.103  a.u.  for  the  NDDO  case. 
It  was  also  found  lhat  the  optimal  C-B-C  bond  angle  is  greatly  overestimated  in 
all  cases  and  the  problem  is  especially  pronounced  when  using  the  LS  resonance 
integrals.  As  discussed  previously  this  problem  can  be  partially  corrected  by 

orbitals.  One  undesirable  side  effect  of  this  scaling  in  the  case  of  borirane  is  that 
the  boron-hydrogen  bond  length  is  greatly  increased  (to  1.36  A). 

Such  an  effect  indicates  that  a single  atom-independent  scale  factor  alone 
will  not  be  sufficient  for  a systematic  improvement  of  geometries  using  the  LS 
resonance  integrals . This  of  course  implies  introducing  atomic  parameters  into  the 
method.  Another  possibility  is  that  the  orbital  exponents  obtained  from  Slater's 
rules  for  INDO  and  used  "as  is"  for  NDDO  and  LS  resonance  integrals  in  this 
work  may  not  be  optimal. 
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Bond/Angle  NDDO  NDDO/M  INDO  INDO/M 


B-C  1.5033  1.5759 

C-C  1.5036  1.4573 


B-H  1.1435  1.0923 

C-H  1.0817  1.0323 

C-B-C  60.014  55.079 

H-C-H  109.42  112.38 

C-B-H  149.99  152.46 

B-C-H  120.72  117.44 


1.5113  1.6139  1.54 

1.5175  1.4728  1.56 

1.1589  1.1047  1.17 

1.0961  1.0459  1.08 

60.269  54.295  60.6 

110.31  113.79  112.8 

149.87  152.85  149.7 

120.47  115.72  120.1 


The  geometries  of  benzene  and  pyridine  were  optimized  with  the  results  as 
shown  in  Tables  16  and  17.  They  are  well  reproduced. 

Table  16:  Calculated  and  observed  equilibrium  bond  lengths  for  benzene6 
Bond/Angle  NDDO  NDDO/M  INDO  INDO/M  Obs.  [93] 
C-C  1.3806  1.3563  1.3942  1.3757  1.397 

C-H  1.0798  1.0312  1.0959  1.0479  1.084 

6.  All  bond  angles  were  constrained  to  the  sp2  hybrid  angle  of  120°  for  the  optimization. 
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Bond/Anglc 

N,-C6 

Cj-Cs 

C4-C5 

CVH, 

Cj-H8 

c2-h, 

CeN.Cj 

C5C4C3 

C,C3C2 

C4C3H8 

NiCjH, 


lengths  (A-B)  and  angles  (A-B-C)  fo 


NDDO  NDDO/M  INDO 


1.3316 

1.3801 

1.3823 

1.0800 

1.0825 

1.0795 

120.31 

119.93 

118.29 

122.35 

115.88 


1.3078 

1.3665 

1.3810 

1.0482 

1.0622 

1.0282 

117.96 

117.76 

118.55 

125.33 

117.70 


1.3479 


1.1001 

1.0955 

118.07 

118.43 

119.44 
120.81 
114.89 


INDO/M 

1.3390 

1.3757 

1.3786 

1.0488 

1.0575 

1.0498 


Obs.  [102] 
1.342 
1.391 
1.398 
1.075 
1.080 
1.085 
116.42 
118.60 
118.36 
121.70 
115.44 


of  the  four  methods  over  all  of  the  molecules  studied  is  presented  in  Table  18. 
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Tabic  18:  Composite  root-mcan-square  and  average  errors  for 
bond  lengths  and  angles  over  all  of  the  molecules  studied. 

Bond  Length  Errors  (Angstroms) 

NDDO  NDDO/M  1NDO  1NDO/M 

RMS  0.065  0.096  0.066  0.097 

AVG  0.011  0.022  0.008  0.019 

Bond  Angle  Errors  (Degrees) 

NDDO  NDDO/M  1NDO  1NDO/M 

RMS  1.657  4.890  1.119  3.316 

AVG -0.198 -0.771 0036 -0.350 


The  success  of  the  INDO/S  method  for  the  prediction  of  electronic  spectra 
is  well  known.  This  being  the  case,  it  would  seem  that  any  other  method  which 
is  to  achieve  widespread  acceptance  must  perform  at  least  as  well  Furthermore, 

there  must  be  enough  of  an  improvement  over  INDO/S  to  justify  this  additional 
cost.  However,  if  there  is  not  a significant  improvement  in  the  results,  the  new 
method  may  still  be  justifiable  on  other  grounds  such  as  a more  sound  theoretical 
foundation  or  a lesser  degree  of  parameterization.  With  the  results  that  follow  we 
hope  to  demonstrate  that  the  methods  we  shall  refer  to  as  NDDO/S  and  NDDO/MS 
not  only  rest  on  a more  sound  theoretical  foundation  but  also  could  be  capable 
of  greater  accuracy  than  INDO/S. 


Configuration  Interaction 


The  calculation  of  electronic  transitions  is  a two-step  process.  The  first 
step  is  an  MO-LCAO  SCF  calculation  yielding  molecular  orbital  coefficients  and 
eignevalues  for  the  Hartee-Fock  ground  state.  These  calculations  require  as  input 

center,  and  the  spin  multiplicity.  Geometries  were  taken  either  from  experiment 
or  the  1NDO  optimized  geometries  when  these  were  deemed  reasonably  close 
to  the  accepted  experimental  geometries.  The  basis  set  is  inherently  valence- 
only,  minimal  basis  with  orbital  exponents  from  Slater’s  rules.  Only  ground  state 
singlet  molecules  were  studied. 

configurations  are  generated  by  single  excitations  from  the  closed-shell  ground 
state  into  the  unoccupied  or  viriual  orbitals.  If  the  ground  state  is  represented 

*0=  (174) 

#5  = (175a) 

= hh7i>»?a7j«4fa7jr} 


(175b) 


(175c) 
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3#J  = I^I?t-0,^c-0AT0«) 

*j=  hh?i-^i?«-^nr7iir)  (175® 

These  are  restricted  open-shell  determinants  which,  in  general,  are  not  eigenfunc- 
tions of  S2,  except  when  all  the  open-shell  electrons  have  the  same  spin.  However, 
by  forming  the  appropriate  linear  combinations  of  those  determinants  which  are 
not  eigenfunctions  of  S2,  new  wavefunctions  may  be  obtained  which  are  eigen- 
functions of  S2.  Such  linear  combinations  are  spin  adapted  configurations.  For 
the  case  of  the  single  cxcilation  described  above,  the  determinants  ' 'I'S  and  34'“ 
are  pure  triplets  as  indicated  by  the  left  superscripts.  Using  the  remaining  two 
determinants,  we  can  foim  the  spin  adapted  configurations 

'*“  = (»;  + *?)  (176) 

= (177) 

where  the  first  is  now  a pure  singlet  and  the  second  is  a pure  triplet.  Considering 
only  singlet-singlet  transitions,  the  elements  of  the  Cl  matrix  are  given  by 

(*o|flr|‘»})  = 0 (178) 

in  accordance  with  Brillouin’s  theorem  and 


= (e.  - e,)6jt6ai  - { ai\bj ) + 2(oi|>6)  (179) 


Spectroscopic  Parameterization 


The  results  which  follow  are  comparisons  among  four  different  methods,  two 
at  the  INDO  level  and  two  at  the  NDDO  level.  The  INDO/S  results  are  well 
known  but  are  included  to  serve  as  benchmark  comparison  for  the  INDO/MS, 
NDDO/S  and  NDDO/MS.  All  four  methods  are  described  briefly  beginning  with 
INDO/S  since  the  others  are  most  succinctly  described  in  terms  of  how  they 
differ  from  INDO/S 
INDO/S 

The  INDO/S  method  of  Zemer  and  Ridley  is  well-known  for  it  accuracy  in  the 
prediction  of  the  spectroscopy  of  large  molecules.  The  essential  elements  of  the 
INDO  method  have  been  described  previously  so  only  the  details  specific  to  the 
spectroscopic  parameterization  need  be  given  and  can  be  summarized  as  follows. 
1.  One-center,  two-electron  integrals  ~iaa: 

7AA  = *3(")  = Ia-Aa  (180) 


2.  Two-center, 


(181) 


(182) 


: naphthalene  i 
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where  the  value  1.2  was  chosen  empirically  to  reproduce  the 
anthracene  spectra. 

3.  Two-center,  one-electron  integrals,  0fAvB\ 

pBA,B  = (183) 

where  S is  the  orbital  overlap,  the  /3°'s  are  chosen  to  give  agreement  with 
experiment  and 

K„  = 1 (184a) 

K,r  = 1 (184b) 

K„  = 1.267  (184c) 

I<„  = 0.585  (184d) 

INDO/MS 

The  INDO/MS  method  is  very  simply  obtained  from  the  INDO/S  method  by 
replacing  the  Wolfsberg-Hclmholtz  formula  for  the  two-center,  one-electron  core 
integrals  with  the  equation  of  motion  formulas  due  to  Lindcrbeig  and  Seamans  as 
described  in  chapter  3.  All  other  tetms  are  as  for  INDO/S. 

NDDO/S 

The  NDDO/S  method  is  a spectroscopic  modification  of  the  NDDO  method 
described  in  earlier  in  this  chapter.  This  modification  consists  of  replacing  a 
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spherically  averaged  piece  of  each  purely  coulombic  two-electron  integral  with 
its  Mataga-Nishimoto  counterpart  to  obtain  a new  coulomb  integral  according  to 
(1‘AI‘aWbvb)'  = [(iiaHaWbvb)  - + 1AB 

(185) 

(J‘a>'aWb*b)'  = [{HaVaWb^b)]* 

where  k is  a scaling  factor  taken  to  be  0.6  where  two-electron  integral  scaling  is 
indicated  and  unity  otherwise.  The  replacement  of  the  averaged  theoretical  integral 
with  its  Mataga-Nishimoto  equivalent  is  done  in  order  to  preserve  the  Mataga- 
Nishimoto  functional  form  for  coulomb  integrals  as  required  for  the  accurate 
calculation  of  tf-ir"  transitions.  The  scaling  of  the  two-electron  integrals  is  based 
on  the  observation  that  the  theoretical  integrals  evaluated  over  the  Slater  basis  are 
approximately  40%  larger  than  the  corresponding  empirical  Slater-Condon  factors. 
NDDO/MS 

This  method  is  the  NDDO  counterpart  of  INDO/MS;  i.e.,  the  Wolfsberg- 
Hclmholtz  resonance  integral  formula  is  replaced  by  the  Linderberg-Seamans 
formulas. 

Spectroscopic  Results 
Formaldehyde 

Unsaturated  molecules  containing  atoms  such  as  nitrogen  or  oxygen,  of  which 
formaldehyde  is  a common  example,  often  undergo  forbidden  n-ir’  transitions. 


: often  the  first  (lowest  energy)  transitic 
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spectrum  they  have  been  of  great  experimental  and  theoretical  interest.  These 
transitions  are  often  poorly  reproduced  by  INDO  methods  which  seems  to  be 
a consequence  of  the  two-electron  scaling  beneficial  to  the  prediction  of  ir-ir‘ 

The  observed  formaldehyde  spectrum  is  shown  in  Table  19  and  the  calculated 
spectra  in  Tables  20  through  22.  From  these  tables  it  is  clear  that  the  INDO/S 
method  does  not  do  a good  job  with  formaldehyde  and  that  the  1NDO/MS  method 
is  not  an  improvement.  The  NDDO/S  and  NDDO/MS  methods  without  scaling 
of  the  two-electron  integrals  seems  to  reproduce  well  the  first  n-ir’  transition  but 
subsequent  transitions  are  rather  dramatically  blue  shifted.  The  NDDO/S  method 
with  scaled  two-electron  integrals  seems  to  represent  the  worst  of  each  extreme 

high.  The  most  interesting  of  all  of  the  calculated  formaldehyde  spectra  is  the 
scaled  two-electron  NDDO/MS  spectrum  which  is  in  the  neighborhood  of  the 
INDO/S  method  for  the  first  n-ir"  transition  but  very  accurate  for  the  next  two. 


Tabic  19:  Observed  formaldehyde  spccirun 


ved  [103] 


28310 

57014 


71541 


a)  The  assignment  of  these  transitions  is  ambiguous 
as  they  seem  to  possess  some  Rydberg  character. 


Table  20:  The  INDO/S  and  INDO/MS  spectrum  of  formaldehyde. 


INDO/S 

INDO/MS 

Type 

tnergy 

Oscillator 

Type 

Energy 

Oscillator 

nV 

25945 

0.00 

n-rr* 

24828 

0.00 

nV 

65250 

0.03 

n-ir* 

50729 

0.01 

rrV 

70979 

0.33 

ir-ir* 

56577 

0.15 

nV 

71584 

0.02 

rr-rr* 

70221 

0.00 

«* 

82896 

0.03 

n-rr* 

70309 

0.04 

103 


Tabic  21:  The  NDDO/S  and  NDDO/MS  (no  scaling  of 


NDDO/S 

Energy 

29123 

85093 

86290 

88090 

105089 


Oscillator 


Typ= 


NDDO/MS 

Energy  Oscillator 
(cm1)  Strength 

28125  0.00 

66589  0.01 

74979  0.42 

82488  0.00 

89815  0.19 


Table  22:  The  NDDO/S  and 
of  the  iwo-clcctron  integrals) 


NDDO/S 

Energy  Oscillator 
(cm1)  Strength 

25367  0.00 

76384  0.02 

77426  0.05 

78371  0.48 

96473  0.04 


24563 

58087 

64916 

80371 

81342 


Oscillator 

Strength 

0.00 


0.01 

0.24 

0.00 

0.12 
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The  ultraviolet  spectrum  of  benzene  is  well  known  and  consists  of  two  primary 
bands  (iw  = 55.900  cm1,  = 48,972  cm1)  and  a secondary  band  at t/(io  = 
38,090  cm'1.  All  three  of  these  bands  correspond  to  rr-tr*  transitions.  The  second 
primary  band  is  a dipole  allowed  transition  (A|9  -*  £|„)  polarized  in  the  x,y 
planes  and  characterized  by  a large  extinction  coefficient  (experimental  oscillator 
strength  of  0.690).  The  first  primary  band  and  secondary  band  correspond  to  dipole 
forbidden  transitions  (A|9  — * B i„  and  .4jy  — » B-ju,  repectively)  and  have  small 
experimental  oscillator  strengths.  The  observed  benzene  spectrum  is  summarized 
in  Table  23, 


Observed  [104] 


Energy  Oscillator 
(cm1)  Strength 

38090  0.001 

48972  0.100 

55900  0.690 

>60000 


The  calculated  INDO/S  spectrum  of  benzene  (36]  (Table  24)  is  in  good 
specifically  for  benzene.  Also  in  Table  24  are  the  INDO/MS  results  which  give 


INDO/S 
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values.  The  most  notable  feature  of  the  INDO/MS  results  is  the  intrusion  of  a 
series  of  ff-tr*  excited  states  between  the  first  and  second  primary  bands  which  is 
a result  of  the  reordering  of  the  degenerate  a mo's  in  the  occupied  space  and  the 
transposition  of  the  highest  it*  mo  with  a o'  mo  (relative  to  INDO/S).  While  this 
could  be  considered  a serious  flaw  of  the  INDO/M  method7,  it  must  be  recalled 


that  the  model  is  parameter  free.  This  considered,  these  results  are  not  too  bad 
and  might  be  brought  more  in  line  with  the  INDO/S  results  by  a simple  scaling 


of  the  calculated  LS  resonance  integrals. 

Table  24:  The  INDO/S  and  INDO/MS  spectrum  of  benzene 


Strength 

0.0001) 


0.0029 

1.0372 

1.0372 

0.0000 


7.  Then 
55,900 1 


for  o-x'  states  below  the 


mat  evid 
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Tables  25  and  26  give  ihe  NDDO/S  and  NDDO/MS  spectra  of  benzene 
without  and  with  scaling  of  the  two-electron  integrals.  The  NDDO/S  results 
arc  clearly  flawed  by  the  poorly  positioned  second  rr-a*  transition.  This  problem 
is  corrected  somewhat  in  the  corresponding  NDDO/MS  spectrum  but  it  is  clear 
upon  comparison  with  the  Table  26  that  the  scaling  of  the  two-electon  integrals 
is  necessary  to  reproduce  the  benzene  spectrum. 

The  NDDO/MS  spectrum  of  Table  26  is  somewhat  remarkable  in  that  it 
represents  a qualitatively  accurate  reproduction  of  the  observed  benzene  spectrum 
using  an  essentially  parameter-free  semi-empirical  model.  Quantitatively,  the 
spectrum  is  blue-shifted  by  approximatley  3,000  cm'1. 
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60481 

60481 

70308 

71536 

71536 

73072 


74863 

74863 

81600 

81600 

81908 


64823 

66722 

68914 


59331 

70947 

70947 

74876 

74876 

75688 


Pyridine 


The  experimental  pyridine  spectrum  is  presented  in  Table  27.  With  the 
exception  of  the  low-lying  n-rr’  transition  due  to  the  presence  of  the  nitrogen 
lone  pair,  it  is  very  similar  to  the  benzene  spectrum.  The  introduction  of  the 

the  Da,  symmetry,  reducing  it  to  C2».  This  perturbation  is  manifested  in  the 
calculated  pyridine  spectrum  as  a splitting  of  the  degenerate  E,„  state  into  two 
non-degenerate  states  of  A|  and  B2  symmetry. 

The  INDO/S  pyridine  spectrum  of  Table  28  gives  generally  good  agreement 
with  the  observed  spectrum.  As  was  the  case  for  benzene  and  the  pyrdine 
spectrum  was  the  basis  for  the  choice  of  so  this  agreement  is  again  not 
surprising.  The  INDO/MS  spectrum  (also  in  Table  28),  on  the  other  hand,  does 
not  correspond  well  with  experiment  and  is  rather  dramatically  red-shifted  relative 
to  the  INDO/S  spectrum.  Furthermore  the  formerly  degenerate  E]„  transitions  of 
benzene  are  now  split  by  a a- ir*  transition  which  seems  to  indicate  that  the 
presence  of  nitrogen  represents  a larger  perturbation  of  the  benzene  system  in  the 
INDO/MS  model  than  one  might  expect. 


27:  The  observed  spectrum  of  pyrdrn 


Observed  [IPS] 

Energy  Oscillator 
(cm  1 1 Strength 

34771  '0.003 

38350  '0.040 

49750  -0.10 

'55000  '1.3 

'56405  diffuse 


Table  28:  The  1NDO/S  and  1NDO/MS  spectrum  of  pyridine 


43917 

49336 

55377 

56067 

57597 

61577 

61939 


25789 

35729 

35730 
45763 
49307 
51876 
52928 


The  results  for  pyridine  using  ihe  NDDO/S  and  NDDO/MS  models  and 
unsealed  two-electron  integrals  are  presented  in  Table  29.  While  the  NDDO/S 
spectrum  gives  the  transitions  in  the  same  relative  order  as  the  INDO/S  model 
the  spectrum  as  a whole  is  blue-shifted  from  6000  to  8000  cm'1  and  the  relative 
positions  of  the  transitions  are  not  well  reproduced.  An  interesting  aspect  of  the 
unsealed  NDDO/MS  spectrum  is  accuracy  with  which  the  first  n-7r*  transition  is 

approximately  4000  cm'1. 

Scaling  the  two-electron  integrals  red-shifts  the  NDDO/MS  spectrum  uni- 
formly by  aprroximately  3000  cm'1  relative  to  the  unsealed  calculation.  The 
NDDO/S  spectrum  is  also  red-shifted  relative  to  its  unsealed  counterpart  but  not 
as  uniformly  with  shifts  varying  from  about  2500  to  4000  cm'1. 


9:  The  NDDO/S  and  NDDO/MS 


electron  integral  scaling)  of  Pyridine. 


42356 

52577 


Oscillator 

Strength 

0.02 

0.07 

0.00 

0.17 


62195 

71702 

72960 

77124 


1.19 

0.05 

0.01 

0.00 


NDDO/MS 

Energy 

34550 

42331 

47138 

54728 

58658 

59614 

67252 

68490 

70132 


Oscillator 

Strength 

0.01 

0.06 


1.07 

1.31 


0.06 


Table  30:  The  NDDO/S  and  NDDO/MS  (with 


llator 

nglh 


02 

00 

01 

05 

19 

00 

00 

02_ 


55859 

56877 

62733 

62931 

66390 


The  observed  furan  spectrum  is  given  in  Table  31.  The  calculated  spectra 
appear  in  Tables  32  through  34.  Many  of  the  trends  already  noted  among  the 
four  methods  are  again  present.  In  this  case  however  the  first  two  observed  bands 
are  best  reproduced  by  the  NDDO/S  and  NDDO/MS  methods,  both  without  two- 
electron  scaling.  These  methods  however,  predict  subsequent  bands,  which  are 
believed  to  possess  some  degree  of  Rydberg  character,  much  too  high.  The 
corresponding  scaled  methods  give  better  overall  agreement  with  the  observed 
spectrum  although  the  first  two  transitions  arc  red-shifted  between  2,000  to  4,000 


Table  31: 


Type 


Observed  [106] 


(cm-1) Strength 


46512  - 
51282 


52219  - 
57143 
58824  - 


sharp  bands 


a)  Possible  Rydberg  . 


n b)  Rydberg  character 


INDO/S  and  INDO/MS  spectrum 


46054 

50171 

52305 

60043 

63694 

64268 

64393 

65074 


0.2456 

0.0055 


43782 

48423 

53656 

55404 

57504 

58300 

60348 

61288 


0.0000 

0.0017 

0.9799 

0.0842 


49663 

54241 

71085 

79855 

83329 

83821 

87084 

91144 

95214 


1.2534 

0.1160 


0.0107 

0.0000 

0.0066 


NDDO/S 


48281 

59658 


67856 

71326 

71796 


47101 

50546 

66027 

69857 

70469 

71470 

73661 

81407 

83546 


CHAPTER  5 
CONCLUSIONS 


The  NDDO  geometries  reported  here  are  only  slightly  better  than  those 
obtained  from  INDO  while  the  INDO/M  and  NDDO/M  geometries  are  not  quite 

INDO  is  a vcty  highly  refined  method  for  which  some  effort  has  been  expended 
in  the  determination  of  the  optimal  set  of  parameters.  Such  is  not  the  case 
with  the  NDDO  results  reported  here  and  1 am  confident  that  a systematic 
reparameterization  of  the  NDDO  method  presented  herein  will  yield  geometries 
consistently  superior  to  those  obtained  from  INDO. 

As  for  INDO/M  and  NDDO/M,  these  are  approximate  methods  in  which  the 
only  parameters,  save  the  orbital  exponents,  are  the  atomic  core  parameters  as 
determined  from  spectroscopy.  That  the  results  are  not  only  reasonable  but  quite 
good  in  some  cases  is  remarkable.  That  they  are  not  as  good  as  the  INDO  and 
NDDO  methods  does  not  necessarily  doom  them  to  the  graveyard  of  unsuccessful 
approximate  methods.  As  is  often  the  case  in  research,  the  stopping  point  of  one 
project  becomes  the  origin  of  several  new  ones  and  there  are  several  aspects  of 
the  equation  of  method  motion  which  warrant  further  investigation. 

The  first  of  the  these  is  the  choice  of  the  functional  fotm  for  on  which 
US 
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Ihe  remaining  core  Hamiltonian  matrix  elements  (resonance  integrals)  depend. 
Linderberg  and  Seamans'  choice  has  already  been  discussed  (see  Eq.  165  on  page 
70).  Other  choices  are  possible.  For  example,  Zemer  and  Parr  (1071  suggested 
a Taylor’s  series  expansion  of  ft,  according  to 

ft.  = + • • • (186) 


Since  they  were  able  to  show  that  ft,  and  B",  can  be  approximated  as  constants 
times  the  appropriate  orbital  exponents,  a one  parameter  (per  atom  pair)  tit 
results  from  which  they  were  able  to  reproduce  force  constants  and  anharmonic 
corrections  for  a wide  assortment  of  diatomics. 

Another  possibility  would  be  to  maintain  the  Linderberg  and  Seamans  choice 

respects  such  a scheme  would  be  preferable  to  that  above  in  that  it  introduces 
no  more  parameters  per  atom  (one  for  an  sfi  basis)  than  the  Wolfsberg-Helmholz 
scheme.  The  rational  being  that  no  more  parameters  are  introduced,  the  formulas 
are  more  soundly  based  in  theory  and,  if  the  results  are  as  good  or  better  than 
those  obtained  from  Wolfsberg-Helmholz,  there  is  no  loss  of  accuracy. 

Also  of  interest,  but  not  explored  in  this  work,  is  a relation  between  the  one- 
center  core  Hamiltonian  matrix  elements  which  is  a further  result  of  the  equation 


method.  This  relation  is 


KA  = <iA{U„  - V„)  (188) 

cannot  be  chosen  independently.  Since,  these  parameters  come  from  atomic 
spectroscopy  and  are  essential  to  a method  which  aspires  to  predict  molecular 

given  the  spectroscopic  C/w' s,  it  is  possible  to  evaluate  Eq.  (188)  for  the 

A preliminary  investigation  of  this  idea  yielded,  for  carbon,  an  exponent  of  1.29 
(compared  to  1 .625),  a value  which  is  not  outside  the  range  of  potential  usefulness. 

Common  to  both  of  the  NDDO  methods  examined  here  is  the  neglect  of 
penetration  approximation  implemented  by  setting  Vab  equal  to  Zb{,basbWa»b)- 
The  observed  tendency  of  the  NDDO  two-center  Fock  matrix  elements  to  be 
greater  in  magnitude  than  their  INDO  counterparts  is  consistent  with  the  overall 
shorter  bond  lengths  and  blue-shifted  spectra.  Since  the  two-center  core  terms 
have  already  been  addressed,  the  logical  next  point  of  attack  could  be  the 
modification  of  Zb^a^bWa^b)  to  Zb(iia»b\«a»b)\I(R)\  where  f{R)  is  some 
function  of  the  distance  between  the 
parameterized  if  necessary. 


This  fu 
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the  use  of  d orbitals  and  the  obvious  extension  to  transition  metals.  The  prospects 
here  are  exciting  and.  for  the  most  part,  unexplored.  Some  work  in  this  area  has 
been  done  using  the  Li  nderbcrg- Seamans  resonance  integrals  [108]  but  only  at  the 
INDO  level  with  a highly  parameterized  spectroscopic  Hamiltonian.  The  effect 
of  the  many  mote  two-electron  integrals  at  the  NDDO  level  on  the  geometries 
and  spectra  of  transition  metal  complexes  is  difficult  to  predict  but  certainly  of 

The  prospect  of  a accurate,  non-empirical,  approximate  method  for  molecular 
electronic  structure  may  seem  like  a glimmer  in  the  eye  of  so  many  hopeful 
theoreticians.  However,  such  a method  would  be  of  enormous  importance  for 
the  application  of  quantum  mechanical  methods  to  large  and  diverse  molecular 
systems.  It  is  my  hope  that  this  work  is  another  step  toward  the  formulation  of 


such  a method 


APPENDIX  A 
THE  IN  DO  METHOD 

Zero  Differentia!  Overlap 

The  differential  overlap  of  atomic  orbitals  on  different 
neglected. 

= S,A„ 

Two-Electron  Integrals 

The  ZDO  approximation  is  again  applied  to  the  two-elec 

WffckfiAp)  = (PAac\’,A>'C)l’AB&CD6p'‘6<r 
However,  among  the  one  and  two-center  integrals,  all  non-coul 
types  are  neglected  while  all  one-center  integrals  are  retained 

it  is  expedient  to  express  the  summation  over  basis  function 


different 


separated.  Therefore, 


with  A.  B,  C,  and  D 


(193) 


B = E EE  ^ac(^aAci^) 
-•£££**»**** 


(194) 


(195) 


(196) 


of  the 
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A = B: 

E EE  E 

+ E E E P-„Xc  (/‘A*cI/M<M)*.4C*»» 

CM  *c  ( 

+ EE  Y.P"^o(m^a\m<’d)6ad6xc 

+EEp«^>''Mw<7m> 

~\  E E E E P’oAc  {m^dm^AdSac  (197) 

-;E  E E p«ic Wc 

'CM  V --  , 

”5  E E E ^^d^a^aWd^a^^x^ao 

- ^EE^-  {PaIaWaPa) 

E EE  J*-Xc*<r<^/lAd/'J4AC7> 

+ E E p’‘»*  [ouXilw*4>  - \(ha*aWaha)\ 


(198) 


2.  Case  11  # v, 
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A = D: 


GHa"a  = 53  53  53  HP«o»c(/MAcKffD)6coWA» 

CfA  Ac  »o 

+ 53  5ZSP»^c0‘AAcl«'AOA>«AcViA„ 

CfA  Ac  

+ 53  y.  ^^SAAo^A^AkA^oJ^AO^^Air 

+ P’a*a  {PA^aWaPa) 

-5  53  53  53  53f»DAc  {iaa^cWdva^ad^ac^^au 

-5  53  5353P<’'>Ac(/‘A-'d<’A‘'A)^C^ciAi. 

~ CfA  Ac  

“513  53  ^PoaXo^A^A^DVA^AD^A,. 

- 5 53 5 3j»a*a  (/aaAaI®a"a) 

(199) 


Gm»a  = 53  53  ^VaAa  [(M^aI*U«a)  - |{/M*/tl<M,M)]  (200) 
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1Z 1Z 1Z  P">Xc  (^/»Acl"a»D>^B*CD 

~ n 1~]  1~!  1^  1^. -P»pAi;{/‘^^dgO''B)^O^BC^^> 


G<MMl  = H H 1Z  1Z  P»D»c</‘/Ackfl‘'fl)^/tD^BC^«4Bi 
= _ 5 5Z  IZ  1Z  P'’D*c(<M‘/cl'70‘'flM/IDiaC<V 
= “ ^ 1Z  1Z 

= ~^lZf»‘D‘-B(^,,fll/'0,'B)^0 

= -^f(M-B{f'^>'fll/<M‘'fl)  (202) 


G|MV0  = -5P|MBB 


(203) 


We  can  now  define  Pec  = 2-PAcAc  10  arr've  at 

1.  Case  n = v\  A = B: 

gimim  = 1Z  PcoOtA^clMA^c) 
cm 

+ 1Z 1Z  P’**«  [(^a^a^a^a)  - 
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2.  Case  /i  ^ r,  A = B: 

= Y.  £ P’^  H»a*a\wa)  - \<j*a>>aY>a»a)]  (205) 


Cmw  = -~P*a*d  {pAvn^Ava) 


Rotational  Invariance 

To  preserve  rotational  invariance  in  the  two-electron  integrals  their  values  are 
assumed  to  be  independent  of  the  nature  of  the  orbitals  Le. 

OmAcI/t-iAc)  = 1AC  {A*c}  (207) 


G|U|U  = 51  PCClAC 


+ £ £ P°**A  [(>,a>>a\ha°a)  - ^a^aWaha)] 


2.  Case  /i  ^ v\  A = B: 

= ££ P"aA. {^a\«a«a)  - \(^aWa»a)]  (209) 
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3.  Case  /x  i/;  A ± B: 


(210) 


The  Core  Hamillonian 

those  of  the  CNDO  method  (see  chapter  2,  section  2). 

Resonance  Integrals 

The  two-center  cote  Hamiltonian  matrix  elements  ate  given  according  to 


and  are  generally  incorporated  into  the  method  as  purely  empirical  parameters. 
The  resulting  equations  are  summarized  at  the  end  of  chapter  2,  section  3. 


= OmI  - -VA-  Vail's)  = /Wa 


(211) 


APPENDIX  B 
THE  NDDO  METHOD 

Zero  Diatomic  Differential  Overlap 


The  differential  overlap  of  atomic  orbitals  on  different 
neglected. 


„„  = 


centers  is 
(212) 


Two-electron  Integrals 

The  ZDO  approximation  is  again  applied  to  the  two-electron  integrals  such 

{paocWbId)  = (213) 

In  doing  so,  for  the  case  of  = v\  A = B,  the  rhs  of  Eqs.  (194)  and  (195) 
vanish  as  does  the  exchange  piece  on  the  rhs  Eq.  (193).  This  leaves 
1.  Case  ii  = v,  A = B: 

+ T,T,p^Ai‘a>‘a\i‘a<’a)  (214) 

- ^ 51 5Z  PuxApa^aWaPa) 


For  the  case  of  /t  ^ w, 


rith  the  obvious  difference  that  /i  and  v are  no  longer  the 


2 Case  /i  f i/;  A = B: 

= Z ZZp»cjcOmvk"c) 

- {/‘A^aWai'a) 


3 Case  fi  ^ v\  A B: 

G-«-b  = ZZZZ  P.nic^^AckflOoJ^flfco 

’ o ' (216) 

“ \ 51  Z Z Z P«DJc(t‘^Ac|oOVfl)^D6flC 

= -jZZ,**»*W»l,'-','»l  <217> 


The  Core  Hamiltonian 


HZh  = M - |v2  -VA-  Vg\uB)  = 


(218) 


APPENDIX  C 

SLATER-CONDON  FACTORS 
Two-Electron  Integrals 


= fl.,(r)y,  ,„(«,*)  (219) 

where  Rnl{r)  is  the  radial  part  of  the  wave  funcdon  and 

Y,Mt)  (220) 

are  the  well-known  spherical  harmonics,  the  general  one-center,  two-electron 
integral  is 


Hi‘>  = J f J ¥>?(lK(2)iW(l)¥.,(2)dT1r,  (221) 

c explicitly 

(*1/0  = J ^K,,.(l)K./.(2)Kn,l,(l)«»,l,(2)r?rJdr1dr.2 
x J J (222) 

xjj  K,;mii(2)y,,mii(2)sin(>2dM<ki 
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The  operator  jL  can  be  expressed  in  terms  of  an  expansion  over  the  spherical 

^ = ^frP,W(c<w  0,)P,W(coS  (223) 

After  substituting  Eq.  (223)  into  Eq.  (222)  and  abbreviating  the  radial  part  of 
the  integral  as 

Rl(i*Jt)  = J fl;,(i(l)K.l.(2)^rB1,,t,(l)«».!1(2)rfr|<lr1dr2  (224) 


{i*m  = £ £ 

x ] JYUf^)yim^)YU l)anfl|<M|^l  (225) 

x / / ^i!mi,(2)K/,m(i(2)Ktm(2)sin  BidBidii 
Now,  integration  over  the  azimuthal  angles  yields  a vanishing  result  unless  die 

m = m,.— m,,  = m,,  - m,.  (226) 

and  reduces  the  summation  over  m to  a single  term.  Furthermore,  the  integrals 
over  the  polar  angle  can  be  expressed  in  terms  of  the  Clcbsch-Gordon,  or  Wigner, 


nple,  Jackson!  1091 


102  or  Slater  [1 10)  page  308. 
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id  Shortley  [111,  112]  as 


.1)/* 

-Ki)T 


/(2la  + l)()2-|m^j)! 

V ('*+KI>! 


'-p 


lmcgrals 


; l,mt,  ij)  (229) 


(230) 
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Fk(nilr,nilj)  = lt(iy,ii)  (231) 

WW)  = f^at(lim,.-,lmi)Ft(nilr,nili)  (232) 

Exchange  Integrals 

{im  = £c*(limil;ljmiJct(lmirJjmil)li‘l(ij;ji)  (233) 

= [ct(l1m(,;f,mlj)J  (234) 

G*(nl/i;nJ/J)  = (235) 

fob')  = (236) 

An  Example 

Consider  the  2p  electrons  of  the  determinant 

|(lOOa)(lOO0)(2OO<*)(2OO/3)(211a)(211/)))  (237) 
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of  the  configuration  (ls22s22pJ).  This  determinate  has  M i = 2 and  Ms  = 0 and 
therefore  must  belong  to  the  term  'D.  The  two-electron  interactions  of  this  single- 
determinant state  can  be  written  (ignoring  the  2s  electrons  since  their  inclusion 
docs  little  to  enhance  the  illustration)  as 

<2p(+|)2p(+i,|2p(+1)2p(+|))  (238) 

i.e.  a single  Coulomb  integral  with  both  electrons  in  a 2p  m,  = 1 spatial  orbital. 
According  to  the  expressions  above  this  integral  simply  evaluates  to 

((211a)(211/?)|(211o)(2n0))  = U)F*(21;  21)  (239) 

Now,  although  the  sum  over  k is  infinite,  it  turns  out  that  there  arc  only  two  values 
for  which  a1  is  nonzero-zero  and  two.  This  leaves 

{(211o)(211/3)|(211a)(211)3)>  = F°(21;21)  + if*(21;21)  (240) 

Note  that  in  this  example,  the  actual  spherical  harmonics  were  assumed  for  the 
angular  part  of  the  function.  In  practice,  quantum  chemists  normally  use  the  linear 

Considering  now  the  exchange  integral  (pip.|psp«>  where 
P«  = 2po 


= ^=(2p(-i>+2p(+1)) 


(241) 


= ^(2p(-i)2po|2po2p(_,)> 

- i(2p(-i)2po|2po2p(tl)) 

- ^(2p(+i)2po|2po2p(_|)) 

- ^{2p(+i(2po|2po2p(+i)) 


(244) 


APPENDIX  D 
MISCELLANEOUS 


Matrix  Elements  of  [r,h] 

Beginning  with  the  definition  of  the  commutator  and  the  assumption  of  a 
complete  basis  set,  we  can  make  use  of  the  resolution  of  the  identity  to  obtain 


- 53omIM*/>)(*dMp») 


Using  the  assumption  2 of  the  development,  namely  that  the  matrix  elements  of 
the  position  operator  are  given  by 


(pa|r|vfl)  = 6ab[6^Ra  + (pair  - /ial"a)] 


(247) 


and  the  one-center  relative  position  operator  matrix  elements  in  an  ,t,  p basis  arc 
(pak  - R-alt'a)  = “alex^A*  + 

e,(^,«x,  + ^,«w)+  (248) 


(pa|r-RaM  = «. 
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where  u*  is  an  atomic  transition  moment  integral.  Then  substituting  Eq.  (247) 
into  Eq.  (246)  and  performing  the  summation  leaves 
(MA|r%B)  = £ 

+ 

(250) 

+ ®si  MumAiam  + 

for  the  first  term  in  the  commutator  and 

= £ {*i.„»DRB(/<A|k|Afl) 

+“fll  e*  (6»R,a6\DtD  + S,„,Si„a) 

(251) 

+ e,(S»..s6x„.+S,n.*x..,) 

+ ei 

for  the  second.  Finally,  summing  over  A in  Eqs.  (250)  and  (251)  gives 
{l‘A\rh\i/B)  = R-aOmIAKb) 

+ “aI®i  (&iuia(za\I‘\‘,b)  + ^iMx<(-S/t|h|"fl)) 

(252) 

+ ®»  (<BaBb(|M|A|i'£)  + 

+ ®.  (^I^M)] 

(n*\hT\vB)  = Rb{ma|A|»'b> 

+ “el®x  (ftwOulty**}  + *Bfl«fl(B4|A|*B» 

(253) 

+ e»  (<B,.«WIAIVi>)  + *BB»s{/M|%fl)) 

+ e.  («.,.b(/mIM*b)  + «,B..MM*fl») 
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(/mIMM  = (R/t  - Ra) WW-'b) 

+ e,  MS,A..{xA\h\uB)  + 6fAX.{,A\h\VB)) 

- + Wl>(/M|fc|3fl))  ) 

+ e,M6w.<{»yt|A|i'B)  + «iu,.(^|fc|«'fl))  (254) 

M^.^aWvb)  + S.„.Mh\sB))  ] 

+ e,  |AM  + Wl^a)) 

«B(6„„,(ax|A|sfl)  + W0(fM|AM)  ] 

Momentum  Rcprcscnlation  of  Linear  Momentum  Matrix  Elements 
In  coordinate  space  a two-center  momentum  operator  matrix  clement  is  written 


P„»  = J *'(r)p*„(r-RK,r 

= J ^(rJjV^r-RJt^r 


(255) 


(256) 


(257) 
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the  general  momentum  matrix  element  becomes 

p„„  = J pXr(p)X*(p)eip'Rd3p  (258) 

The  plane  wave  e'p'n  can  be  expanded  [113]  as 

e""R  = E E c>'<"  j'Kpfl)  1T(*.  « (259) 

and  in  this  case  it  is  sufficient  to  specialize  the  expansion  to 

e'pR  = ;'i(pR)  fl(coafi)  (260) 

$ being  the  angle  between  p and  R.  Using  the  orthonormality  of  the  Legendre 
polynomials  we  find  that 

“I  jl(pR)  = — j—  J e"Kco§9Pt{cos0)dcosO  (261) 

which  can  be  compared  with 

il(pR)  = 27?  | eifRcom9  Pi(cos  9)d  cos  0 (262) 

to  find 

eipRto*9  = (21  + l).1  j,(pX)  P,(cos0)  (263) 

Now,  using  this  expansion  in  the  expression  for  the  gradient  matrix  element  gives 

Ppv  = J pXp(p)Xp[p)  elpR  d3,, 

= jpXp(p)Xr(p)'£(2l+l)i,jl(pR)Pl(cOS9)d3p  (264) 

which  is  useful  for  deriving  relationships  among  the  momentum  matrix  elements. 
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Identities 

v"  = 1r(RV")  (265) 

Beginning  with  the  evaluation  of  the  left-hand  side 
= Xa.(p)p^xbm^^p 

= £ JxA.(p)  ^sin  0 cos » r+  sin B sin  (.]+ M*B.(p) cosfle,>R<Pp 
(266) 

Terms  1 and  2 vanish  as  can  be  seen  by  symmetry  arguments  considered  earlier 
or  by  explicit  integration  over  the  azimuthal  angle.  Considering  now  the  third 

P»  = V„e,8 

V"  = i / / / $**«**.«  “*S  9 e'^Vdp  sin  «<»</<)  (267) 

and  integration  over  the  azimuthal  angle  gives 

V,’  = kJ  J WXa.(p)Xb.{p)“*'‘  6 e"‘ndpsinfl,M  (268) 

Now,  for  the  left-hand-side 
v„  = ^ J XA,(p)^XBw(p)  sin  ft cos  (S  e'>*R<Pp 

= ^ J XA.(p)^  fsinfleo^r  + sinflsin^J  + co^jj  A'fl.(p)  sin 0 cos  <J  c""R#p 

2 3 (269) 
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p„  = V„e,  (270) 

V"=^/ / / (P)^fl. (P) sin2 « c‘pn cos2 « pepsin flrffld,! 

= \J  f \p\Xa.(p)Xb.(p)  sin’  9 c,p  “dpsin  9d9  (271) 

Since  it  is  our  intention  to  differentiate  V„  with  respect  to  R,  it  is  expedient  to 
do  so  before  expanding  the  plane  wave  to  arrive  at  the  final  expression.  However, 
we  will  first  rewrite  Eq.  (271)  as 

v”  = i // IpI^p)*".  (p)  si“2  6 r'pKdpS\a  6d6 

= \ j J lpl^.(p)7fB.(p)(l  - co s29)  e'-X^'dp  d(cosfl) 

= '-lJ  (1 -cos'fl)  «"■*“*»  d(cos  9)  (272) 

I = J \p\Xa.(p)XbM  ^ (273) 

Now.  by  the  product  rule 

^tK-)-«^v.  + v.  ««) 
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= x2)ipRx  c'fR‘  dx  (275) 

where  x has  replaced  cos  8 for  brevity.  Now.  substituting  Eqs.  (272)  and  (275) 
into  (274)  gives 

= j/  |x(l  - xJ)  ipR  ci,K*  dx 

+ ^lj(l-  x1)  ci,Rz  dx  (276) 

After  integration  by  parts  the  first  integral  on  the  right-hand  side  of  Eq.  (276) 
becomes 

'-I  J (3xs  - 1)  ei,R’  dx  (277) 

Finally,  combining  the  intermediate  results  gives 

= ‘-i  j (3x2  - 1)  e**  dx  + Li  j (1  - xs)  e'>&  dx 
= ‘-I  J x2  e,pfll  dx  (278) 
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which  after  replacing  x with  cos#  ar 


he  definition  of  / le 


^(HV„)  = j \P\Xa,{p)Xb,{p)  cos2  # 'i,R“ 

= '0  } J |p|Jf-t.(p)Xfl,(p)eosJ9e1> 
Finally,  comparison  of  the  above  with  Eq.  (268)  gives 


(280) 
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